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USAF  Contract  No.  AF  18(600 )-90.  The  contract  is  administered  by 
the  Mapping  and  Charting  Branch  of  the  Photographic  Reconnaissance 
Laboratory,  Weapons  Components  Division,  Wright  Air  Development 
Center,  Wright-Patterson  Air  Force  Base,  Ohio.  Mr.  A.S.  Rosing  is 
Project  Engineer  on  the  project  applicable  to  the  subject  of  this 
report. 

Research  and  Development  Order  Nos.  R-6S3-44,  ’’Research  in 
Photo grammetry  and  Geodesy  for  Aeronautical  Charting,"  and  R-683-5^, 
"Aeronautical  Charting  Systems,"  are  applicable  to  this  report. 

This  report  was  originally  initiated  at  the  Ohio  State  University 
Research  Foundation  as  OSURF  Technical  Paper  No.  159* 
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ABSTRACT 


In  the  present  Technical  Paper  are  given  the  outlines  of  Intrinsic 
Geodesy,  a  theory  "by  which  the  gravity  field  of  the  Earth  may  he  analyzed 
by  using  only  the  results  of  actual  measurements,  and  without  any  ad¬ 
ditional  assumption  on  the  structure  of  the  field  itself.  Any  geometrical 
or  mechanical  element  that  would  not  have  a  physical  reality  is  therefore 
excluded. 

After  having  specified  the  aims  of  the  problem  and  the  means  by 
which  it  may  be  solved  (Introduction  and  Section  I),  short  summaries  on 
the  elements  of  vector  and  tensor  calculus  are  given  (Section  II  and 
HI). 

By  the  help  of  this  calculus,  the  study  of  the  gravity  field  of 
the  Earth  is  performed,  using  coordinates  innate  in  the  field  itself 
(Section  IV),  and  several  applications  to  practical  problems  are  shown. 

Somigliana's  field  is  finally  suggested  as  the  most  appropriate 
for  generalizing  the  usual  ellipsoidal  Geodesy  in  a  three  dimensional 
scheme  (Section  V). 
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the  Air  Force  of  the  findings  or  conclusions  contained  therein.  It 
is  published  only  for  the  exchange  and  stimulation  of  ideas, 

FOR  THE  COMMANDING  GENERAL: 


Chief,  Pho/to  Reconnaissance  Lab 
Weapons  Components  Division 


WADC  TR  52-149 


iii 


TABLE  OE  CONTENTS 


SECTION  PAGE 

INTRODUCTION  _ _  _  1 

I  PRELIMINARY  _  5 

II  VECTOR  CALCULUS  (ELEMENTS  AND  NOTATIONS)  _  12 

III  TENSOR  CALCULUS  (ELEMENTS  AND  NOTATIONS)  _  27 

IV  INTRINSIC  GEODESY  _  52 

V  APPLICATION  OF  INTRINSIC  GEODESY  _  75 

REFERENCES  84 


WADC  TR  52-149 


iv 


INTRINSIC  GEODESI 


INTRODUCTION 

As  for  arty  other  Science,  the  development  of  Geodesy  has  been 
greatly  influenced  by  the  historical  process;  in  our  Science  it  led 
us  from  the  primitive  idea  of  a  flat  Earth  to  the  actual  concept  on 
the  shape  and  size  of  our  planet,  and  its  internal  structure. 

In  early  times  it  was  discovered  that  the  Earth  could  be  approxi¬ 
mated  by  a  sphere,  and  already  the  Greeks  and  the  Arabs  tried  to  deter¬ 
mine  its  radius.  Only  much  later  the  theoretical  works  of  Galileo, 
Newton,  Jhygens,  and  Clairaut  led  to  the  ellipsoidal  conception  of  the 
shape  of  the  Earth  that  was  definitely  proved  ty  the  celebrated  ex¬ 
peditions  of  the  Acadeny  of  France  in  Lapland  and  in  Peru  as  early  as 
the  second  quarter  of  the  18th  century. 

As  geodetical  work  progressed  and  methods  and  instruments  were 
greatly  improved,  it  was  soon  discovered  that  even  the  ellipsoid  could 
only  be  regarded  as  a  further  approximation  for  the  shape  of  the  Earth; 
and  the  idea  of  the  Geoid  (Listing,  1873)  became  essential. 

Notwithstanding,  and  since  the  times  of  the  expeditions  of  the 
Acadeny  of  France,  the  ellipsoid  was  still  the  dominating  feature  in 
Geodesy;  not  only  it  was  and  is  still  universally  used  in  practical 
Geodesy  and  surveying;  but  it  also  affects  most  of  theoretical  research. 

It  may  be  pointed  out  at  once  that  the  function  of  the  ellipsoid 
in  Geodesy  may  be  regarded  as  two-fold.  Sometimes  it  is  accepted  as  a 
good  approximation  to  the  Geoid,  and  therefore  confused  with  it.  That 
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happens  mostly  in  problems  of  practical  Geodesy,  such  as  triangulation. 
Sometimes  it  is  simply  regarded  as  a  standard  surface  of  reference,  with 
respect  to  which  the  "anomalies*1  or  the  "deflections"  of  the  Geoid  are 
referred.  But  several  times  its  use  is  promiscuous,  and  it  is  not  easy 
to  discern  what  in  the  results  is  due  to  the  actual  structure  of  the 
Earth's  gravity  field,  and  what  is  affected  by  the  ellipsoidal  assumption. 

It  seems,  therefore,  of  some  interest  to  establish  a  theory  leading 
to  the  study  of  geodetical  and  geophysical  problems,  abstracting  from 
any  hypothesis  on  the  shape  of  any  level  surface  in  the  gravity  field 
of  the  Earth,  and  leaning,  therefore,  only  on  the  data  of  actual  measure¬ 
ments. 

Moreover,  since  the  time  Geodesy  did  no  longer  confine  to  the  study 
of  geometrical  problems  on  the  shape  and  the  size  of  one  particular  level 
Surface,  and  started  studying  the  dynamical  problems  of  gravity,  its 
concern  became  more  ample,  and  its  aim  more  precise.  Geodesy  is  the  Science 
devoted  to  the  study  of  Earth's  gravity  field.  The  study  of  the  geoid 
became  a  particular  problem  in  this  general  definition. 

Our  study  must,  therefore*  not  be  confined  to  one  particular  level 
surface*  but  must  extend  to  the  third  dimension  in  such  a  way  as  to 
link  in  an  overall  picture  both  the  geometrical  and  the  dynamical  aspects 
of  geodetic  problems. 

The  difficulty  to  overcome  at  first  lies  in  the  fact  that  Operative 
Geodesy  must  always  refer  to  coordinates  suitable  for  carrying  out  actual 
measurements;  if  independence  from  any  exterior,  arbitrary  system  of 
reference  (like  the  ellipsoid,  or  cartesian  baricentric  axes,  etc.)  is 
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claimed:,  then  the  coordinates  to  b@  used  must  be  innate  in  the  field 
itself.  We  will  therefore  call  such  coordinates  intrinsic,  and  our 
study  Intrinsic  Geodesy. 

Also,  in  such  study  we  will  only  consider  entities  having  a 
physical  reality  and  quantities  susceptible  to  undergo  actual  direct  or 
indirect  measurements  which  Operative  Geodesy  is  already  able  to  perform, 
or  would  be  able  to  perform  if  appropriate  theoretically  conceivable 
instruments  should  be  available.  A  standard  surface  of  reference  (like 
the  ellipsoid)  or  a  standard  field  of  reference  (like  Somigliana' s 
field)  should  be  sometimes  used  for  testing  the  actual  gravitational 
field  of  the  Earthy  but  th®  particular  assumption  of  such  standards 
would  never  affect  the  general  results. 


We  already  remarked  that  the  ellipsoid  may  be  correctly  regarded 
as  a  standard  reference  surface  in  two-dimensional  Geodesy  for  referring 

the  discrepancies  between  it  and  the  Geoid.  The  most  natural  extension 
of  such  concept  in  three  dimensions  is  given  by  Pizzetti-Somigliana's 
field  of  gravity,  which  was  accepted  in  1930  by  the  Stockholm  meeting 
of  the  International  Association  of  Geodesy  on  the  suggestion  of 
Dr.  Lambert,  Prof.  Cassinis,  and  Prof.  Somigliana  himself.  Somigliana' s 
field  is  fully  determined,  on  account  of  the  celebrated  Stoke' s  theorem, 
by  giving  the  ellipsoid  as  one  of  its  level  surfaces  (usually 
Hayford's  International  Ellipsoid),  by  the  well-known  angular  speed  of 
rotation  of  the  Earth,  and  by  the  equatorial  value  of  gravity  on  the 
ellipsoid  itself  (usually  g0  *  978.OI4.9  gal).  On  Somigliana' s  field  is 
based  the  international  formula  for  normal  gravity,  the  coefficients 
of  which  have  been  determined  by  Prof.  Silva  and  Cassinis. 


WADC  TR  52-149 


3 


The  principles  of  Intrinsic  Geodesy  may  be  very  easily  particularised 
for  Somigliana' s  field,  and  the  results  used  as  a  standard  of  reference 
or  as  a  good  first  approximation  in  whatever  geodetic  problem  the  third 
dimension  is  involved;  and  that,  in  a  perfectly  similar  way  to  that  in 
which  the  ellipsoid  is  used  in  two-dimensional  Geodesy. 

Modern  techniques  tend  indeed  each  day  more  to  extend  geodetic  ex¬ 
ploration  in  space,  as  is  shown, e.g.,  by  the  impressive  fact  that  in 
classical  Geodesy  only  measurements  along  geodesics  and  lines  on  a 
surface  (Geoid  or  Ellipsoid)  were  considered,  and  that  at  present  times 
optical  and  electronic  devices  allow, instead, the  measurement  of  distances 
along  lines  in  space,  which  may  be  very  closely  confused  with  geodesics 
in  space,  i.e.,  straight  lines. 
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INTRINSIC  GEODESY 


SECTION  I 

PRELIMmKI 


1.1  Ob .1  act  of  Geodesy 

We  win  first  point  out  the  object  of  Geodesy,  as  it  is  referred 
to  in  the  present  lectures  %  it  is  the  study  of  the  gravity  field  of  the 
Earth. 

As  it  is  very  well  known,  the  field  of  gravity  (  •  gravitation  + 
centrifugal  forces),  in  the  same  way  as  all  conservative  (non-dissipativ@) 
fields.  Bay  be  thoroughly  described  by  means  of  a  single  apace  function, 
the  potential  function  as  introduced  first  by  Laplace,  in  all  points 
exterior  to  Hi®  volume  occupied  by  attracting  masses.  The  object  of 
Geodesy  may  be  also  formulated,  therefore,  as  the  study  of  the  potential 
function,  or  also  of  the  shape  and  size  of  its  level  surface®  (aquipotantial 
surfaces),  and  of  their  orthogonal  trajectories  (lines  of  force). 

It  is  also  very  well  known  that  the  study  of  a  potential  field 
involves  geometrical  and  mechanical  aspects £  for  instance,  the  problem 
of  determining  the  size  and  shape  of  each  of  two  near  equipotential 
surfaces  is  of  geometrical  nature,  but  that  of  specifying  th®  valu®  of 
the  corresponding  potential  difference  (th®  work  to  be  done  by  trans¬ 
ferring  the  unit  mass  fro®  the  one  surface  to  the  other)  is  a  mechanical 
one,  which  involves  the  concept  of  "force”  (gravity). 

We  may  add  that  in  this  scheme  the  aim  of  Topography  (Surveying) 
is  merely  the  study  of  the  shape  of  the  actual  physical  surface  of  the 
Earth,  referred  to  the  above  gravity  field,  or  in  a  more  simple  way,  to 
a  particular  level  surface  of  the  same. 
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L.2  Tiro  Aspects  of  the  Gravity  .Field  of  the  Earth 

The  study  of  a  potential  field  may  be  considered  from  two 
completely  different  points  of  view.  The  first  (also  historically)  is 
based  on  Newton’s  law  of  gravitation  and  on  elementary  laws  of  centrifugal 

forces^  i.e.,  it  is  based  on  the  consideration  of  the  "cause**  of  gravity 

/ 

itself.  It  admits,  therefore,  an  interaction  between  attracting  masses 
and  is  called,  therefore,  by  Hermann  Weyl  \J  the  "Fernwirkungsgesetz** 
(principle  of  action  in  distance).  Its  fundamental  equations  are,  as 
it  is  easily  seen. 


grad  W  -  vw 


(1-1) 


where  W  is  the  potential  of  gravity  in  a  point  P  of  the  space,  due  to 
the  attracting  masses  of  density/*.  distributed  in  the  volume  (V),  and 
to  a  rotation  of  angular  velocity tu;  f  is  Newton's  constant  of  gravita¬ 
tion,  r  the  distance  between  P  and  any  point  whatever  of  (V),  R  the 
distance  of  P  from  the  axis  of  rotation.  Moreover  g  is  the  vector  of 
gravity  (the  gradient  of  the  potential  i).  (Here  and  in  the  following, 
over lined  letters  like  g  will  stand  for  vectors.) 

We  may  also  refer  to  the  above  principle  as  an  integral 
principle,  owing  to  the  fact  that  the  knowledge  of  the  whole  distribution 
of  attracting  masses  is  required  for  computing  the  integral  in  (1-1) 

It  is  the  principle  followed  in  the  classical  works  of  Clairaut, 
Helmert,  Bruns,  etc. 

The  second  principle  may  be  referred  to,  on  the  other  hand,  as  a 
differential  principle,  no  longer  involving  the  knowledge  of  the  distri¬ 
bution  of  attracting  masses.  Weyl  refers  to  this  principle  as  to  the 


Hermann  Weyl,  Raum,  Zeit»  Matcrie  -  Vorlesungen  liber  a] 
Relativitdtstheoris.  (Berlin,  1921).  (Also  in  French.) 
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"Wahewirkungsgeset*"  (principle  of  action  in  proximity).  This  denomina¬ 
tion  is  due  to  the  fact  that  only  the  knowledge  of  the  field's  structure 
in  the  immediate  neighborhood  of  a  given  point  is  wanted. 

In  other  words,  the  gravity  field  may  be  studied  here  as  an  "effect", 
independently  from  its  "cause"  (the  distribution  of  masses) *  and  it  is 
well  known  that  the  gravity  vector  g  must  satisfy  the  following  two 
fundamental  (differential)  equations: 

rot  g  ■  v  x  g  *  0  ;  div  g  «  v  •  g  m  2 ax 2  -  h^f/*  .  (1-2) 

The  first  equation  affirms  that  the  field  is  conservative,  and 
the  second  that  a  rotation  is  superposed,  and  that  in  general  in  the 
point  P  where  g  is  considered,  masses  of  density/*  are  present  (the 
last  term  disappears  in  vacuo,  and  to  the  greatest  approximation  in 
free  air). 

The  differential  principle  has  been  followed  ty  Stokes  primarily, 
and  his  theorems  are  of  utmost  importance  in  Geodesy*  and  later  on  ty 
Pizzetti,  Somigliana,  and  many  other  prominent  geodesists  and  physicists. 

The  outstanding  importance  of  the  second  principle  lies  in  the 
fact  that  it  is  completely  independent  of  the  distribution  of  attracting 
masses,  and  it  may  therefore  be  applied  with  great  advantage  in  Ge¬ 
odesy,  where  the  distribution  of  density  in  the  interior  of  the  Earth 
is  unknown. 

1.3  Geodetic  Position  of  the  Problem 

We  have  so  far  outlined  the  problem  of  gravity  from  the  stand¬ 
point  of  mathematical  physics.  Our  task  is  now  to  examine  it  from  the 
point  of  view  of  Geodesy,  i.e.,  of  experimental  science,  and  to  specify 
suitable  methods  for  the  practical  exploration  and  numerical  definition 
of  the  field  itself. 
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On  account  of  its  independence  of  the  knowledge  of  the  Earth's 
interior,  in  the  following  we  will  confine  ourselves  to  considering 
only  the  second  principle  stated  above. 

Our  problem  may  be  simply  formulated  as  the  problem  of  studying, 
in  a  given  domain,  the  vector  £  of  gravity,  and  other  vectors  connected 
with  it  in  a  very  simple  way.  Such  vectors  are,  for  instance,  the  unit 
vector  n  of  g  (g  -  -g  71),  where  g  (  a  scalar)  is  the  intensity  of 
gravity,  or  "gravity";  the  derivatives 

dg  d*g  dn  dan 

as  9  as2  *  *  *  *  i  dS  *  dSa  *  *  *  * 

of  the  same  vectors  in  given  directions,  the  vector 

Q  -  grad  g  (1-3) 

(the  gradient  of  gravity),  etc.,  and  the  only  difficulty  is  to  find  a 
way  of  defining  numerically  the  above  absolute  (not  connected  with  a 
particular  reference  system)  entities. 

A  vector  itself  is  not  a  measurable  entity.  Only  its  components 

(or  projections)  with  respect  to  a  given  system  of  three  cartesian 

axes  are  numerically  definable  and  therefore  apt  to  be  measured  by 

means  of  suitable  instruments  and  methods,  arid  to  be  used  in  numerical 

(or  algebraic)  computations.  are  therefore  forced  to  chooses  in 

Geodesy,  reference  systems  of  three  axes  (orthogonal  or  not*  unit  or 

not,  but  not  coplanar)  in  respect  to  which  to  consider  the  components 

(or  the  projections)  of  our  vectors,  and  it  is  a  question  of  fundamental 

importance  for  the  practical  use  in  Geodesy  to  have  in  mind  that:  ‘ 

U)  the  axes  of  reference  must  have  an  immediate  physical  reality; 

(2)  no  additional  hypothesis  must  be  made  in  connection  with  the 
gravity  field  itself; 
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(3)  the  components  (or  projections)  of  all  vectors  to  be 
considered  may  be  easily  measured-* 

The  three  conditions  above  exclude  therefore,  for  instance,  an 
absolute  cartesian  system  of  reference  for  the  whole  Earth,  which  would 
not  have  a  physical  reality  and  would  make  it  therefore  impossible  to 
measure  coordinates  and  components-# 

So  far  as  the  possibility  of  having  a  single  system  of  reference 
for  the  whole  Earth,  satisfying  the  above  conditions,  is  obviously  to 
be  excluded,  we  are  obliged  to  choose  suitable  local  reference  systems, 
one  for  each  point  of  the  field  in  the  region  to  be  considered. 

As  a  first  solution,  reference  systems  of  this  kind  would  be 
given  by  orthogonal  unit  systems  of  axes  in  any  point,  one  axis  being 
directed  towards  the  Zenith,  a  second  towards  the  North,  and  a  third 
towards  the  East.  Such  systems  are  often  used  in  local  geodetic  surveys 
of  limited  extension,  and  we  will  refer  to  them  as  local  astronomical 
systems.  Each  aystem  satisfies  obviously  the  three  fundamental  condi¬ 
tions  stated  above  and  has,  therefore, an  intrinsic  definition* 

But  we  must  always  have  in  mind  that  the  extension  to  a  wider 
domain,  as  is  often  considered  in  Geodesy,  requires  the  different  systems 
of  reference  to  be  connected  each  to  the  other  in  such  a  way  as  to  make 
it  possible  and  easy  to  transfer  all  mathematical  properties  from  on# 
system  to  the  next,  or  also  to  canpare  the  components  (or  projections) 
of  the  same  entity  referred  to  different  points  in  space* 

We  will  see  that  the  local  astronomical  systems  are  not  convenient 
for  such  purposes,  in  spite  of  their  very  simple  definition,  and  this 
is  due  to  the  fact  that  it  is  not  easy  to  connect  them  each  to  another. 

We  will  see,  on  the  contrary,  that  another  system  is  much  more  suitable 
for  the  same  purpose,  though  its  definition  is  a  little  more  complicated 
than  the  foregoing* 
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1.U  Object  or  Intrinsic  Geodesy 

The  task  of  our  study.  Intrinsic  Geodesy,  may  be:  we  must 
define  in  any  point  of  our  gravity  field  a  system  of  three  non-coplanar 
axes,  each  satisfying  the  condition*  stated  in  1.3  ani  give  a  role 
as  single  as  possible  for  connecting  each  system  to  the  next  one. 

If  we  are  able  to  do  this,  we  will  also  be  able  to  develop  every 
geodetic  research  on  an  actual  basis  without  referring  to  any  arbitraxy 
hypothesis  (like  that  of  the  spheroid  or  the  ellipsoid),  nor  to  any 
System  of  reference  not  accessible  to  our  measurements* 

An  easy  way  to  reach  this  result,  as  we  will  see  in  more  detail, 
is  to  establish  in  the  whole  gravity  field  an  intrinsic  system  of 
curvilinear  coordinates,  i.e.,  coordinates  whose  definition  is  implicit 
in  the  field  itself  and  do  not  need  therefore  any  additional  assumption. 
We  will  see  that  the  most  suitable  coordinates  to  be  used  for  this 
purpose  are  those  defined  by  the  following  coordinate  surfaces: 

(1)  the  family  of  equipotential  surfaces  themselves;  we  at¬ 
tribute  to  each  surface  as  parametric  coordinate  the  corresponding 
value  of  the  potential; 

(2)  the  family  of  the  surfaces  connecting  all  points  in  space 
having  the  same  astronomical  latitude  if  i  i3  at  the  same  time  the 
parametric  coordinate  connected  with  this  family; 

(3)  the  family  of  the  surfaces  connecting  all  points  in  space 
having  the  same  astronomical  longitude  X;  X  is  at  the  same  time  the 
parametric  coordinate  connected  with  this  family. 

The  coordinate  lines  are  therefore: 

(1)  the  isozenithal  lines  (intersection  of  4>  and  X  surfaces), 
connecting  in  space  all  points  having  the  same  Zenith  (same  $  and  X); 
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(2)  the  geodetic  meridians  (intersection  of  W  and  X  surfaces)} 

(3)  the  geodetic  parallels  (intersection  of  W  and  $  surfaces). 

All  entities  so  defined  have  a  physical  reality;  the  equipotential 

surfaces  can  be  materialized  by  a  level;  astronomical  latitudes  and 
longitudes  may  be  also  easily  measured;  the  difference  of  potential  W 
is,  with  abstraction  of  an  insignificant  numerical  factor,  the  dynamic 
difference  of  height. 

Attention  is  recalled  to  the  fact  that  geodetic  meridians  and 
parallels  should  not  be  confounded  with  North  and  East  tines,  i. e.,  the 
lines  whose  tangents  are  directed  in  each  point  towards  the  North  or 
East,  and  also  that  geodetic  parallels  and  meridians  do  not  cut  in 
general  at  right  angles. 

As  it  may  be  easily  seen,  the  system  of  coordinates  so  defined 
is  intrinsic;  in  fact,  we  did  not  use  any  additional  hypothesis,  or 
any  exterior  element  other  than  the  axis  of  the  World  to  which  astro¬ 
nomical  latitudes  and  longitudes  are  referred.  Moreover,  the  direction 
in  space  of  the  axis  has  a  physical  reality  and  may  be  used  for  ex¬ 
perimental  measures. 

The  coordinate  lines  so  defined  give  us  now  the  possibility  of 
fixing  the  position  of  any  point  in  space  by  means  of  its  three  para¬ 
metric  coordinates  $,  X,  and  W.  And  moreover,  the  coordinate  lines 
defined  above  enable  us  to  draw  at  each  point  a  system  of  three  funda¬ 
mental  vectors  tangent  to  the  lines  themselves.  Furthermore,  as  we  will 
see  in  more  detail  afterwards,  the  same  coordinate  lines  will  allow  us 
to  connect  the  systems  of  fundamental  vectors  in  neighboring  points  and 
enable  us,  therefore,  to  compare  vector  entities  at  points  far  apart. 
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SECTION  II 


VECTOR  CALCULI®  (ELEMENTS  AND  NOTATIONS) 

2.1  Preliminary 

In  the  following  we  shall  consider  absolute  entities,  like 
points,  vectors,  and  holographies,  which  all  are  independent  from  any 
system  of  reference;  and, on  the  other  hand,  relative  entities,  like 
coordinates,  components,  and  projections,  which  all  depend  not  only 
upon  the  absolute  entity  they  represent,  but  also  upon  a  reference 
system  of  coordinates. 

The  absolute  entities  are  studied  by  the  Vector  Calculus,  and 
the  relative  ones  by  Cartesian  Geometry  in  the  elementary  case  of 
rectilinear  reference  systems  of  coordinates,  and  by  Tensor  Calculus 
in  the  case  of  general  curvilinear  coordinates.  (See  reference  U ,  5, 

6,  11.) 

The  difference  of  two  points  P2  -  P1  *  3  defines  a  vector;  we 
may  write,  therefore,  also  P2  »  Px  +  3.  From  this  point  of  view,  a 
may  be  also  regarded  as  an  operator  leading  from  Px  to  P2. 

A  vector  is  completely  defined  by  a  direction  (of  the  straight 
line  joining  Pj  and  P2),  a  sense  on  it  (from  Px  to  P^),  and  a  magnitude, 
or  modulus,  or  length  (of  the  segment  P,  P„).  A  vector  may  be  reoresented 
by  an  arrow,  ah  vectors  specified  by  the  same,  three  qualities  stated 
above  are  said  to  be  equipollent.  The  point  from  which  the  vector  starts 
is  immaterial. 

We  will  always  use  an  overlined  letter  for  indicating  vectors. 

The  same  not  overlined  letter  will  represent  the  modulus  (length)  of 
the  vector  (a  scalar). 


A  unit  vector  is  sometimes  called  a  versor;  so, for  the 

_  5 

instance,  u  ■  —  is  a  unit  vector;  we  have  therefore  "S  ■  a  ET.  A 

vector  may  also  be  defined  by  its  modulus  and  its  unit  vector. 

Let  us  have  twa  infinitesimally  near  points  P  and  P* ,  and  put 

P*  -  P  »  dP;  dP  is  therefore  an  infinitesimal  vector,  and  we  may  say 

that  the  differential  of  a  point  is  a  vector.  If  ds  be  the  length 

PP’  then  £  .  —  is  a  unit  vector.  If  P  and  P*  are  two  infinitesi- 
ds 

mally  near  points  of  a  curve,  then  F  is  the  tangent  unit  vector  to 
the  curve  in  P. 

We  assume  that  the  operations  of  sum  and  subtraction,  and 
multiplication  by  a  scalar  of  vectors  are  well  known. 

Scalar  product  of  two  vectors  "S  and  F  is  indicated  by 
a.b  (Italian  notation,  a  x  F)  (read  5  dot  F,  or  also  W  scalar  F) 
and  gives  a  scalar  (a  pure  number).  If  9  be  the  angle  between  the 
(positive  directions  of  the)  two  vectors,  then  we  have 

a.b  »  a  b  cos  0 

Scalar  or  dot  product  is  distributive  and  commutative. 

Condition  for  orthogonality  of  two  vectors:  a.b  »  0 
Vector  or  cross  product  of  two  vectors  is  indicated  by 
ax  d  (Italian  notation,  a  A  F)  (read  a  cross  F,  or  also  a  vector  F). 
For  full  definition  of  vector  product  a  preliminary  orientation  of  space 
is  necessary.  If  we  are  given  tnree  (non-coplanar)  vectors  vx,  v2,  v3, 

(in  the  given  order),  such  orientation  is  fixed  by  the  rule  that  an 
observer  situated  along  one  of  the  three  vectors  would  see  the  other 
two  (considered  in  cyclic  order),  the  one  at  his  left*  the  second  at 
his  right*  following  one  of  the  two  schemes: 


WADC  TR  52-149 


13 


In  the  first  case  we  will  say  that  the  orientation  of  our 
space  is  dextrorsum*  or  clockwise,  whereas  in  the  second  case  it  is 
sinistrorsm,  or  anticlockwise.  We  will  assume  for  our  space  clock¬ 
wise  orientation. 

Vector  or  cross  product  of  two  vectors  *  and  E  is  a  new 
vector  C  the  direction  of  which  is  perpendicular  to  both  the 
directions  of  H  and  E.  Its  sense  is  such  that  (a,E,E)  be  a  posi¬ 
tive  trihedron  (dextrorsum  in  our  case),  and  its  modulus  ab  sin  0. 

Vector  product  is  distributive,  but  not  commutative: 

axE«-Exa 

Condition  of  parallelism  of  two  vectors:  a  x  E  «  0. 

Triple  product,  or  mixed  product,  or  box  product  of  three 
vectors  a,  E,  c  is  the  scalar 

V  -  [a  E  c]  -  axE.c-cxa.E-Exc.a« 

-  a.Exc»c.axE-E.cxa* 

*  -b*a,c  « . 

We  give  here  some  useful  elementary  formulae: 


(II-D 
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(a  x  F)  x  c-a.cE-B.cI, 

S  x  (S  x  £)  •I.fC-i.n, 

(a  x  b)  x  c  +  (b  x  7)  x  a  +  (c  x  a)  x  S  -  0, 

(a  x  E)  ,  (f  x  3)  «  a  .  c  a. 3 

_  _  _  _  » 

b .  c  b  .  d 


u .  a 

u  .  "5 

u .  c 

[a  1) 

c]  [u  V  w]  ■ 

v  .  a 

v.  B 

V  ,  c 

% 

w.s 

w.B 

w.  c 

(II-2) 

(a  x 

b)  x  (c  x  d)  < 

-  (a  c 

d)b  - 

[b  c 

d]a  * 

■  [a  b  d]c  -  [a  b  c]  d  , 

a  [E  c  d]  -  b  [c  d  a]  +  c  [d  a  b]  -d[aEc]  -  0  , 

.  td  5  ?]  ;  ♦  [a  d  c]  b  *  t? i ?] «  a;5-]t0). 

[a  Hj  c] 


Condition  of  coplanarity  of  three  vectors:  [a  b  cj  -  0. 

Rules  for  differentiation  of  products  are  the  same  as  for  ordinaiy 


produc  ts : 

d(a  ,"B)  *  a.dB  +  5.  da, 

d(a  x  b)  ■  a  x  db  +  da  x  b  ,  (H”3) 

d[a  b  c]  ■  da  xb.  c  ♦  a  x  db,  c  +  ax  b.dc. 

2.2  Cartesian  Components:  Differential  Formulae 

If  we  are  given  a  positive  system  of  three  orthogonal  unit 

vectors  (i^  Ia,  I3)>  we  W  always  write 

a  *  ajij^  +  a2±a  ♦  a3i3  , 

g^t  SLj,,  a3  are  the  components  (and  also,  in  this  case,  the  projections) 
of  a  on  the  three  axes:  a^.  •  a  .  iy  . 


(II-U) 
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W«  hare: 


b  -  a  b  ♦  a  b  ♦  a  b  :  ax  b  = 
*a  3  3  3  * 


ii  ia  is 

*1  *3  s 

b,  b* 

3  2  3 


( H-5) 

*1  *2  *3 1 

[I  B  c]  -  b*  ba  b3 

I  Cl  C2  c3  j 

Positional  vector  of  a  point  P  referred  to  an  origin  0  and  the 
syste*  (i3>  ia,  i3): 

P-°  -  x *  x2ia  +  x3i3  ,  ( II-6) 

(x1,  Xj,  x3)  are  the  cartesian  coordinates  of  P, 

A  function  W(x^,  x#,  x3)  -  W(P)  is  a  function  of  position. 

We  call  gradient  of  ¥  and  write  grad  W  or  also  vf 
(read  nabla  1)  the  following  vector* 

grad  W  »  vl 


aw  -  aw  -  aw  - 


»X3 


(II-7) 


’•grad*  ia  an  operator  between  scalars  and  vectors. 

It  is  immediately  seen  that  grad  W  .  dP  -  dW;  it  follows  that 
grad  W  is  always  perpendicular  to  the  surfaces  W(P)  -const. 

If  we  are  given  a  vector  W  -  S(P)  function  of  position  (the 
components  of  which  are  functions  of  P),  we  shall  call  curl  or 
rotation  of  M  and  write  rot  ¥  or  also  v  x  u  the  following  vector 

7  „  au  _ 


ixx 

tilt; 

fl**1 

♦ 

X 

*x( 

aXa  »X3  j 

ii 

i. 

9 

a  a 
*x2  ax3 

*1 

u2  u3 

16 

'  9x, 


axji 


ax- 


-3) 
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"pot"  is  an  operator  between  vectors  and  vectors. 

It  follows  immediately  that 

rot  grad  W  «  Vxv  W  -  0  (II-9) 


A  field  of  vectors  u  ( P)  such  that  rot  u  ■  0  identically 
is  called  irrotational  or  conservative.  In  such  case,  a  scalar 
function  W(P)  exists,  such  that  u( P) “grad  W(P).  W(P)  is  called 
scalar  potential,  or  potential  of  the  field. 

If  we  are  given  a  field  of  vectors  u(P),  we  shall  call 
divergence  of  u  arri  write  div  u  or  also  v  .  u  the  following  scalar: 


div  u 


V.  u 


<?U  3u, 

U  +  +  _ 3 

9xJ  ^X2  dx3  '' 


( II- 10) 


"div"  is  an  operator  between  vectors  and  3calars. 


We  have  identically  div  rot  v  *  0;  and,  conversely,  if 
div  u  ■  0,  it  follows  u  *  rot  v. 

A  field  such  that  div  U  »  0  identically  is  called  solenoidal* 
it  follows  that  every  solenoidal  field  is  the  rotation  of  a  new  field 
of  vectors  (vector  potential). 

Combining  the  two  operators  div  and  grad  we  put 

div  grad  W  *  &2w  -  V.VW  »V2W,  (II-ll) 

the  new  operator  A2  (between  scalars  and  scalars)  is  called  the 
Laplacian  operator. 

A  given  field  u(P)  may  alw^rs  be  written  as  the  sum  of  an 
irrotational  and  a  aolerioidal  field.  We  have 

u  -  grad  W  +  rot  v  ,  (11-12) 

and  also  it  may  be  written  as 

u  ■  grad  W  ♦  m  grad  n  ,  (11-13) 

where  W,  m,  n,  and  V  are  suitable  scalar  or  vector  functions  of  P. 
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It  is  very  important  to  point  out  at  once  that  although  the 
foregoing  differential  operators  (grad,  div,  rot)  have  been  defined 
for  the  sake  of  simplicity  by  the  help  of  cartesian  coordinates,  they 
are  absolute  operators,  as  we  can  see  from  the  following  vector 
formulae  which  may  define  as  well  the  operators  themselves: 
grad  W  .  dP  -  dW  , 

rot  u.dPxSP*du.$P-&Tt.dP,  ( II- lU) 

(div  u)dP  x  jP  .  dP  -  «Px  ap  .  d!I+  9P  x  dP  .  Str  ♦  dPx  SP  .  , 

where  d,  6,  d  are  differentials  in  three  non  parallel  nor  coplanar 
directions. 

Following  formulae  may  be  derived  immediately  from  the  foregoing: 

gfad(V  ♦  W)  ■  grad  V  +  grad  W,  rot  (u  +  v)  -  rot  u  +  rot  v, 
grad  niff  «  m  grad  V  *  W  grad  m  , 

rot  m  tl  *  m  rot  u  ♦  grad  m  x  *[ ,  ( 11-15) 

div  mTT  »  a  div  u  ♦  grad  m  .  u, 

div  u  x  v  «  w .  rot  u  -  u  .  rot  v  . 

2.3  Homographies 

A  linear  operator  (transformation)  a  between  vectors  and  vectors, 
such  that 

a  a  -  u  , 

a('ff  +  b)  *  «  a  +  a  b  ,  (distributive  property)  ( 11-16) 

a(m  S)  ■  maM  ,  (commutative  property) 

will  be  called  a  homograpby. 

Each  homograpby  transforms  parallel  vectors  into  parallel 
vectors,  and  vectors  parallel  to  a  plane  into  vectors  parallel  to  a 
plane. 
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Moreover,  each  homography  may  be  fully  determined  by  the  corre¬ 
sponding  vectors  of  three  (independent-non  parallel  to  a  plane) 
given  vectors. 

Proper  homography  is  a  homography  transforming  vectors  non 
parallel  to  a  plane  into  vectors  non  parallel  to  a  plane. 

Singular  or  degenerative  homography  is  a  homography  transforming 
a  group  whatever  of  independent  vectors  into  other  vectors  parallel 
to  a  plane,  or  even  parallel  to  a  direction. 

Some  special  homographies  are: 

(1)  homothety  (proper  homography) 

a  •  m  (real  number),  ax  *  m  x  (x  arbitrary).  ( 11-17) 

(2)  axial  homography  (singular  homography  transforming  a  group 
whatever  of  independent  vectors  into  a  group  of  vectors  parallel 
to  a  plane) 

a  ■  i«  ,  ax-uxx  (II-18) 

(3)  dyad  (of  Gibbs)  (singular  homography  transforming  each  vector 
into  a  vector  parallel  to  a  direction) 

a  ■  H  (u ,  v)  ,  ox«(u.x)v  (11-19) 

(U)  dilatation  (proper  homography) 

Homography  satisfying  to  the  condition 

5c  .  a  y  -  y  .  ax  .  (11-20) 

A  homography  may  he  indicated  also  by  the  following  notation, 
which  defines  it  fully: 

/u  v  w\  a  &  -  u 

a»\__/i.e.,  _  _ 

'a  b  c  '  ab  »  v  (11*21) 

a  c  ■  w  . 
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In  the  case  of  a  dilatation,  there  may  always  be  found  three 
unit  vectors  (principal  directions)  ii>  *2*  ia  such  that 

/mix  *I3  *I3\ 

o-  {  _  j  (11-22) 

V  ii  **  **  /* 

The  three  invariants  of  an  homography  are  Ijtt,  Iaa,  I^a 
(first,  second  and  third  invariant)  and  are  defined  by  the 
following  formulae; 

Ixa  .uxv.w»vxw.au+wxu.av  +  u  x  v  .  a  w  } 

Iaa  .u  x  v  .  w  -  (ay)  x  (aw)  .  u  +  (aw)  x  (au)  .  v  +  (au)  x  (a  v)  .  w  , 

I3a  „  u  x  v  .  w  «  (a*u)  x  (ay)  .  (aw)  .  (11-23) 


The  vector  Va  of  a  homography  o  is  defined  by 

2V  a„  ii  x  v  «v.  alT  -  u  .  av  (u,  v,  w  arbitrary  vectors).  (11-21;) 

The  dHatatim  Da  of  a  is  defined  by 

Da  -  a  -  Va  x  .  (11-25) 

The  conjugate  Ka  of  a  is  defined  by 

Ka  -  Da  -  Vax  .  (11-26) 

The  cyclic  Ca  of  a  is  defined  by 

ca  *  i^a  -  a  .  (11-27) 

Theorem  1;  Da  is  always  a  dilatation. 

Theorem  2;  A  homography  may  always  be  written  uniquely  as 

the  sum  of  a  dilatation  and  an  axial  homography,  as  follows: 

a-Da  +  Vax  t  Ka  »  Da  -  Vax  .  (11-28) 

We  have 

fVa-.o  , 

a  is  a  dilatation  if 

Ka  »  a  ; 
fDa  ■  0  t 

a  is  axial  if 

Ka  ■  -a  . 
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9ooe  important  formulae: 

Ix(a  *f>)  -  Ixa  ♦  Ia/* 

lima 

»  alxa 

V  (a  ♦  p  )  •  va  ♦  lp  , 

Vma 

-  m  Va 

D  (a  ♦/>)  «  Da  ♦  D p  , 

Dna 

-  nDa 

K  (a  ♦  p  )  »  Ka  ♦  K  p  , 

Kma 

■  mXa 

C  (a  ♦/»)*•  Ca  +  C/9  , 

Cma 

■  mCa 

(11-29) 


Ix(u  -  o  , 

V(u  *  v  )  -  u  , 

I1H(u,  v)  -  u.v  , 

2VH(u,  v)  -  5  x  v  , 

KH(u,  v)  ■  H(v,  u) , 


I*®  -  3»a  ,  Ijt.  -  ** 

Efcn  *■  »  ,  K*  -a 

Ia(uxv  )-  u3  ,  I3(uxv).  0 

D(u  x  v  )  ■  o  ,  K(u  x  v  )-  -u 
I2H(u,v)  -  0  ,  IaH(u,T>-  0 
2DH(us  v)  »  H(xjjY)  ♦  H(v,u) 


(a-scalar) 


axialu 


(n-30) 


dyads 


(11-32) 


Commutation  theory  of  Jacobi: 

x.ay  -  y.Kai  ( 11-31) 

(a,  x,  and  y  arbitrary  homography  and  vectors). 

Ii  (/»«)  -  ♦  Xaa  ♦  Ia/,  -  Ia(a  ♦  />)  , 

T  /  \  T  a  T  (II“32) 

ia  (p®)  -  i3a  W  • 

2.U  Spatial  Derivatives  of  Fee  tors 

Jw 

*e  will  define  the  operator  —  and  call  it  spatial  deri- 

dP 

vative  of  the  vector  u  -  u(P),  as  the  homography  which  applied  to  a 
displacement  if  t  gives  us  the  corresponding  increment  STT  undergone 
by  TX: 

%  .  (11-33) 
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to  a  given  unit  vector  a. 


If  we  apply  the  operator  — 

du  _  _ 

—  a  gives  us  the  derivative  of  u  in  the  direction  a.  For  this 

dP 

du  _ 

reason  "jp  may  be  called  the  spatial  derivative  of  u. 


We  have  the  following  absolute  expressions  for  rot  u  and 


div  u: 


—  du 
rot  u  -  2V^  , 


dS 

div  u  -  Ix~  . 


(11-34) 


Furthermore,  the  following  important  formulae  may  be  noticed: 

_  —  du  -  d7  _ 
grad  u  .  v  -  K”  v  +  K~  u  . 

If  ix,  i3,  i3  be  a  constant  orthogonal  system  of  vectors, 
u  ■  u^ij^  +  uaia  +  u3i3  , 


(11-35) 


r°t  u  -  I  grad  *  ir  f 


div  u  -  i  grad  ur  .  ir  . 


(11-36) 


Some  theorems  follow: 

Theorem  1:  If  everywhere  in  a  given  field  u  .  rot  u  ■  0,  then 
three  scalar  functions  of  position  m,  n,  and  p  exist,  such  that 

u  «  m  grad  n  , 
rot  ( fi  u)  •  0  » 

(p  «  1/m)  . 


(II-37) 


Theorem  2:  If  everywhere  in  a  given  field  div  u  -  0,  then 
two  scalars  m  and  n  exist  such  that 

u  ■  grad  m  *  grad  n 


(11-38) 
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Theorem  3:  For  an  arbitrary  vector  field  u,  scalars  m  and 


n  exist  such  that 

rot  u  *  grad  m  x  grad  n  . 

2.5  Cartesian  Components  of  Homographies 

We  may  consider  the  cartesian  components  of  a  homography. 

If  (i^)  are  the  orthogonal  unit  vectors  of  our  reference  system  and 
a  the  homography  considered,  we  may  put 

aix  ■  axi  ^-1  +  aaa  ^2  +  ^3  * 

aTa  -  a21  T1  +  a23  T2  +  a23  T3  > 

^3  “  a3^  *1  +  a32  T2  +  a33  T3  * 

The  matrix  ||  aw  ||  gives  us  the  components  of  a  homography 

and  defines  it  completely.  In  fact,  if  we  put 

u  -  i.  u_  Ir  ,  we  will  have 
r»l  x 

afll  *  u3  alx  +  UjjttTg  +  u3al3  - 

“  («x*ix  +  U2a2X  +  ♦ 

♦(UxSa  +u2a22  +  V»)*a  + 

Vi.  +Uaa23  ♦  V.3>T3  * 

We  have,  therefore: 


(II-39) 


(II-UO) 


( Il-Ul) 


(II-U2) 
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!•  hftTtf  furthermore: 


•» 

*13  *lj 

a  « 

*33  *33 

1 

*3» 

***  *33 

Hi 

*31  *31 

Ka  « 

*13 

Si  *33 

t 

*13 

*33  *33 

Si 

V2(Ha  ♦  *») 

l/2(al3  +  a31) 

Da  ® 

l/2( 

*21  +  *12) 

*33 

V2(*a3  +  *33) 

l/2( 

*31  +  S3) 

V2<«„  *  *aa> 

*33 

2Vo  s 

(*33 

“  *32)^1 

(*3i  "  S3)Ta  + 

^*13  ”  *21  )*3 

0 

l/2( aia  -  a3l) 

V2(S3  “  *3l) 

Vax  a 

1/20 

Si  “  S2) 

0 

V2(aa3  -  *52) 

1/Zh 

i3l  _  S3) 

V2(a3a  -  a23) 

0 

I 

0 

+*3 

■*3 

lx  s  I 

“*3 

0 

hs 

1 

I 

+*3 

-*1 

0 

H(u,v)  = 


ss 

ss 

V. 

Vi 

ss 

ST, 

v» 

v» 

Vl 
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9 


( 11-1*3) 


lift  a  I1Ka 


ix 


+  s. 


28 


+  & 


33 


a32  a23 

a33  a3X 

ail  ®12 

•f 

+ 

*33  a33 

aX3  ^1 

Si  &23 

an 

*12 

a!3 

Vs 

a22 

&23 

a3l 

a32 

S3 

a  -  2- 

a 

H(i  , 

i  ) 

r,s 

rs 

r* 

8 

■  IaKa 


f 


(decomposition  of  a  homograph/  into  dyads) 


Furthermore: 


»UX 

du3 

t 

1 

5xT 

•1 

'•SIS 

/du  3u 

f  aTSa  5x? 

du  \ 

351  . 

OX2 

dU2 

rr- 

0X2 

du3 

3x£ 

* 

I 

l1*  r- 

*»/ 

dua 

du3 

si; 

du 

—  x  ■ 

dP 

♦  X, 

—  du 

> 

du  _ 

K—  x  “ 
dP 

-  du  — 

x.  r~  ix 
dx1 

—  du  - 

+  x  •  dXg  i* 

+  X. 

an  _ 
dx3  i3 

» 

rot  u  * 

—  du 

*x*  y*T 

♦4*1^ 

+  I3 

x  23L 

*  dX3 

» 

div  u  » 

T  du 

Z1  •  " — 
axl 

♦  I 

+  H  •  dx2 

+  I3 

du 

’  aX3  “ 

dUj^ 

^T  + 

dU2  dU3 

5xJ  + 

Ag  W  ■  div  grad  W  « 

T  d  gyad  *  . 

1  dP 

a*w 

♦  S- 

*v 

a*w 

(II -Mi) 
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aaw 

dx^9xa  (Hessian) 


(Jacobian) 
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SECTION  III 


TENSOR  CALCULUS  (ELEMENTS  AND  NOTATIONS) 


3.1  General  Rectilinear  Reference  Systems 

We  have  studied  so  far,  absolute  entities  like  vectors  and 


homographies,  and  also  absolute  operators  like  grad,  rot,  div. 


,  etc.,  either  without  using  any  reference  system  of  axes  or  coordi¬ 
nates,  or  only  using  cartesian  orthogonal  unit  systems. 

As  we  will  see,  it  is  instead  often  very  important  to  refer  the 


same  absolute  entities  and  operators  to  general  cartesian  systems 
(i.e.,  not  unit,  and  not  orthogonal  systems),  as  happens  in  the 
infinitely  small  domain  around  a  point  using  general  curvilinear  coordi¬ 
nates.  (See  reference  3,  4,  6,  7,  8,  9,  10,  11,  12.) 

As  a  first  step,  We  shall  confine  our  study  to  rectilinear 
cartesian  coordinates,  neither  orthogonal  nor  isometric,  as  they  may 
be  specified  by  a  point  0  (the  origin)  and  three  (not  coplanar  and  not 

unit)  vectors  vx,  va,  va,  our  fundamental  (or  base) vectors. 

If  we  are  given  an  arbitrary  vector  U  and  we  want  to  define 


u  by  means  of  our  reference  system,  we  always  may  express  it  as  a 

linear  combination  of  the  fundamental  vectors,  as  follows: 

u  ■  ■*  u*T_  ♦  u*t_  ■  S  u®V 

*  3  8  8 


(III- 


(lower  case  Roman  indices  have  here  and  in  the  following  the  range  from 
1  to  3)*  or  also  we  may  define  it  by  means  of  the  three  scalar  products 

ux  -  TX.  Vx  f  yxz  -  u.  7a  }  u3  -  u .  ?3  .  (III- 

Both  pysteas  of  numbers  ur  “  (u1,  ua,  u3),  or  also 

«a,  u3)  allows  us,  with  the  help  of  our  reference  system 
of  vectors,  to  reconstruct  u$  ur  or  also  ur  are  therefore  called  the 
components  of  u  in  our  system. 
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More  precisely,  ur  are  called  the  contravariant  components 

of  u  ,  Uy  the  COY  a  riant  components.  Contravariant  and  covariant 

components  are  in  general  different.  They  coincide  only  in  the 

elementary  case  of  cartesian  orthogonal  unit  reference  systems. 

The  reason  for  the  denominations  covariant  and  contravariant 

given  to  the  two  kinds  of  components  is  due  to  the  fact  tnat,  by 

transforming  the  original  fundamental  vectors  vr  into  other  new- 

fundamental  vectors  v*  by  means  of  a  linear  transformation 

r 

+  ^2  +  -5a?Vs  >  (III-W 

the  two  kinds  of  components  of  the  same  vector  do  not  transform  by  the 
same  rule,  as  we  will  see  at  once. 

(Summation  Convention).  We  will  use  in  the  following  the  con¬ 
vention  to  dispense  with  the  sign  2  >  all  times  an  index  is  repeated 
in  a  monome  twice,  once  as  a  superscript,  and  once  as  a  subscript. 

Thus,  for  instance,  we  will  write  the  following  formulae s 

u  -  u^g  ;  ? J.  -  a®vs  .  (III-5) 

The  letter  used  for  the  summation  index  (dummy  index)  is 
immaterial  (like  s  in  the  foregoing  formulae),  and  apy  letter  may 
be  substituted  at  will  for  it.  An  index  which  is  not  repeated 
(like  r  in  the  foregoing),  is  called  a  free  index. 

Let  us  now  have,  in  the  first  system  of  reference 
u  -  u^v  ,  u^  -  u  .  vr  ,  and  u  «  u'  rfcj,  ,  u^,  -  u  .  vj.  in  the  second; 

then 

u  -  USVS  ■  U,rvJ.  *  u’  r  (III-6) 

and  therefore 

us  *  u* r  a®  (law  of  contravariance)  ,  (III-7) 
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and  furthermore 


«  u.  v^.  ®  a*  u6  vg  .  a^«s  (law  of  covariance),  (III-8) 

As  we  see,  u^,  transforms  like  the  fundamental  vectors 
v_  themselves  (they  are  therefore  cogradient  to  the  same,  and  the  law 
of  transformation  is  co variant);  the  ur  transform  conversely  by  a 
different,  contragradient  law  with  respect  to  the  fundamental  vectors, 
i.e.  they  are  contr&variant  to  the  same  vectors. 

In  a  quite  similar  way  we  may  consider  the  components  of  a 
homography  a,  and  put 

a  "  *s  (covariant  components) .  (III-9) 


It  may  easily  be  seen  that  these  components  transform  by  the 
following  generalized  law  of  covariance: 


P  <1 

ttr  as  *pq  • 


(in-  io) 
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(we  identify  the  indices  differing  by  multiples  of  3);  and  reciprocally 


-  m£1L±jl2L 

Vr  [v*  v2  v3  ] 


(111-14) 


Moreover 


[v1  V2  V3]  = 

[vi  vz  v3] 


(ui-15) 


By  the  help  of  the  reciprocal  vectors  we  also  may  write 


u»u/ 

r 


r  _  _r 
u  *  u  .  v 


(111-16) 


In  a  similar  way  we  may  put 


rs  _r  s 
a  *■  av  ,  v 


(contravariant  components  of  a) 


r.  _r  _ 

a  «  crv  .  v 

.*  s 


.s  _s 

a  »  ov  .  v 
r.  r 


(mixed  components  of  a)  . 


By  a  transformation  of  the  fundamental  vectors  v  ,  the 

r  ' 

reciprocal  vectors  transform  contravariantly. 


(111-17) 


We  may  put 


s  s_p 
v*  -  A  v 

P 


(III-18) 


and  have 


—  _  s 

v*  .  v' 
r 


~s  t  .  s_  p  ts_p  t  s 

S  *  a  A  v,  ,  /  *  a  A  S  *  a  A  , 
r  rpt  rpt  r  t 


(HI-19) 


and  finally 


t  s  _  s 
°  A  a  S 
r  t  r 


(II 1-20) 
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is  the  inverse  of  the 


matrix 


This  formula  shows  that  ||  A  [j 

t 

||  ||  .  Our  formulae  read  now: 


_  s_ 

v*  -  a  v  , 

V 

r_ 

*  A  v* 

r  r  s 

s 

s 

to 

*1 

(covariance)  , 

u'  *  a  u  , 

r  r  r  * 

u 

s 

-  W 

( III-21) 

r  r_s 

V  *  A  v  , 

-s 

V 

.  «v 

( c  ontravarianc  e)  . 


,r  .r  s 
u*  =  A  v  , 

s 


3.3  Tensors;  the  Fundamental  Metric  Tensor 

r 

3.3.1  All  systems  (single,  double,...)  of  members  like  u^,  u  , 

rS  V* 

a  a  a*«...  having  the  meaning  specified  above,  connected 
rs  *  .s’ 

to  the  fundamental  system  of  vectors  (or  its  reciprocal),  and 
transforming  accordingly  to  the  above  rules  of  covariance  or 
contravariance,  are  called  tensors  (of  order  one,  of  order  two,  etc.). 

If  we  are  given  two  vectors 

u  «  u  vr  «  uSv  ,  w  ■  w  v1*  ■  wSy  ,  (III-22) 

r  s  r  s 


we  may  form  their  scalar  product  in  several  forms: 

_  rsrs  scr  rrsSr 

U.W  -uwg  -uwg  »  U  W  o  “  u  w  =  u  W  O  *■  u  W 

rs  rsrsr  sr  r  y 


having  put 


g 

rs 


_r  -s 
v  .v 


sr 

g 


( 111-23) 


It  may  easily  be  seen  that  g^g  and  grs  are,  respectively,  a 
covariant  and  a  contravariant  double  tensor.  We  call  them  the  funda¬ 
mental  metric  tensor  in  covariant  and  in  contravariant  form. 
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We  also  have 


_  _r _ r  s_  _  _r  s_r 

u  -  u  v  ■  u.v  v  ■  u  v  .r  ▼  -  g  u  v  » 

r  r  s  r  sr 


and  therefore 


u  -  g  u 
r  sr 


In  a  similar  way  we  may  find 


u 


s  r 
g  au  ,  u 
rs 


rs 

g  u  ,  a 
a  rs 


„  Dpq 
g  g  a  , 
rp  sq 


( III-210 


(111-25) 


(III-26) 


The  fundamental  metric  tensor  allows  us  to  compute  the  scalar 
product  of  two  vectors,  and  therefore  to  establish  a  linear  and  angular 
metric;  moreover,  it  allows  us  to  "change  the  variance"  of  the 
components  of  a  given  tensor,  i. e. ,  to  deduce  the  covariant  components 
from  the  contravariant  ones,  and  conversely. 

We  also  have  the  formulae: 


t  “or 


r  s 

3  r 

*  v 
r  r 


_s 

g  v  » 

rs 


rs_ 

g  vr> 


(III-27) 


ys  rs 

i.e.,g  »  o  ,g  ,  are  respectively  the  covariant,  mixed,  and  con- 
rs  i* 

travariant  components  of  the  fundamental  system  of  vectors  and  its 
reciprocal. 


3.1*  Operations  on  Tensors.  Invariants 

Addition_of_ tensors*.  The  addition  of  the  corresponding  components 
of  two  tensors  of  the  same  kind  and  order  gives  another  tensor  of  that 
order  and  kind;  thus  for  instance: 


*  +  b 
v  r 


r  r 
a  ♦  b 
st  st 
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Product  of  tensors.  The  product  of  the  components  of  (two) 
tensors  of  any  kind  or  order  gives  another  tensor  the  order  of  which 
is  the  sum  of  the  orders  of  the  two  given  tensors;  thus  for  instance; 


s 

a  b  *  c  ,  ab 
r  s  rs  r 


s 

c 

r 


pq  n  pqn 
a  b  *  c 
rst  m  rstm 


Contraction  of_ tensors.  If  we  make  an  index  of  contravariance 

and  an  index  of  covariance  the  same  in  a  given  mixed  tensor,  so  that 

it  becomes  a  dummy  index,  and  we  must  sum  for  it  from  1  to  3,  the 

result  is  a  new  tensor,  the  order  of  which  is  reduced  by  two;  thus; 

rst  st  rst  t  rst  t 
a  ■  c  ,a  «  d  ,  a  *  e  etc. 
rm  m  rs  sr 

Invariants.  If  we  contract  completely  a  tensor  of  even  order 
with  as  many  indices  of  covariance  and  contravariance,  the  result  would 
be  an  invariant  (a  scalar  function  of  position,  the  value  of  which 
does  not  depend  upon  the  fundamental  system  of  axes);  thus; 


r  rs 

a  *  b  ,  a  *  c 

r  rs 


rs 

a 

sr 


d  . 


Contracted  jnultiplication.  *e  may  multiply  two  tensors  of  any  order 
and  contract  at  the  same  time  with  respect  to  a  couple  (or  more  than  a 
couple)  of  indices.  The  result  is  a  tensor,  or  an  invariant,  thus: 
a  br  ■  c  (■  a. b,  scalar  product  of  a  and  b;  invariant). 


a^a  ■  (a)2  (square  of  the  length  of  a;  invariant), 

r  s  s  r. 

a  ‘b  «  c  (homography  cc  of  mixed  components  a  applied  to  vector 
•  s  .s 

—  r—  _  s 

b  *  b  v  *  vector  c  »  c  v  1 

r  s'  ’ 

P  r 

a  b  -  c  etc. 
rs  pt  st 


33 
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Quotient  law  of  tensors.  If  the  result  of  the  contracted 

multiplication  of  an  arbitrary  tensor  by  a  given  system  gives  us  a 

tensor  or  an  invariant*  then  the  given  system  is  a  tensor. 

Thus,  for  instance,  if  we  are  given  a  system  of  the  third  order 

s 

A(r,s,t)  and  we  take  an  arbitrary  tensor  b  ,  and  we  know  that 

r 

A(r,s,t)b®  =*  c1  , 

where  c^  is  a  tensor,  then  A(r,s,t)  is  also  a  tensor  and  must  therefore 
be  written  as 

A(r,s,t)  *  art 
s 

The  foregoing  rule  allows  us  to  recognize  the  tensor  character 
of  a  given  system.  It  is  very  important  to  stress  that  the  tensor 

g 

used  as  factor  (in  our  instance  b  )  must  be  completely  arbitrary. 

r 

A  particular  case  is  to  ruse  as  test  tensor  the  product  of  tensor* 
of  the  first  order;  thus,  for  instance,  if 

X*  s 

A(r,s, t)x  y  z  -  c  (an  invariant) 

ti 

A(r,s.t)  is  a  tensor  and  must  be  written  as  a  .We  may  also  say 
*  *  rs 

that  the  coefficients  of  a  multilinear  invariant  form  in  the 

r  s 

variables  x  ,  y  ,  z are  the  components  of  a  tensor. 

We  may  also  observe  that  the  contraction  of  a  tensor  of  the 
first  order  and  the  fundamental  (or  reciprocal)  vectors  give  us  a 
vector,  i.e.,  an  absolute  entity;  thus  ^v,  ■  TI . 

3.5  General  Curvilinear  Coordinates 

In  the  foregoing  we  have  merely  used  cartesian  rectilinear 
coordinates,  i.e.,  the  same  system  of  reference  for  the  whole  space 
considered.  Thus,  for  instance,  if  we  were  given  a  vector  a  the 
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components  of  which  would  be  a*  (or  a^)  at  a  given  point  P  of  space, 
and  we  transport  by  equipollence  the  vector  a  to  another  point  0,  the 
components  of  ¥  would  remain  unaltered.  The  same  may  be  stated  for 
tensors  of  any  order  and  kind. 

In  mary  instances  it  is,  however,  not  possible  or  convenient  to 
use  cartesian  rectilinear  coordinates.  Let  us  confine  our  attention  for 
the  sake  of  simplicity  to  a  planej  it  is  always  possible  to  use  on  it 
cartesian  rectilinear  systems  of  reference,  and  the  foregoing  theories 
may  be  very  easily  applied  with  some  slight  modifications  (summation 
from  1  to  2  instead  of  from  1  to  3,  obvious  modification  as  to  the 
reciprocal  vectors).  But  this  would  no  longer  be  possible  if  we 
consider  a  curved  surface  (not  developable)  on  which  no  cartesian 
system  of  coordinates  is  possible. 

But  even  in  three  dimensional  space  it  may  happen  that  we  were 
led  to  use  curvilinear  coordinates  (for  instance,  spherical  coordinates, 
cylindrical  coordinates,  ellipsoidal  coordinates,  etc.),  as  may  be 
suggested  by  the  particular  problem  we  are  dealing  with.  In  more 
advanced  geometrical  theories  we  may  even  consider  curved  three 
dimensional  spaces  in  which  no  cartesian  system  of  reference  would  be 
possible. 

In  such  cases  the  foregoing  simple  theories  cannot  be  applied 
and  we  must  establish  other  rules  for  the  definition  of  the  components 
of  our  absolute  entities,  like  vectors,  homographies,  etc.  in  any 
point  of  space,  and  for  operation  on  them. 

Let  us  consider,  for  the  sake  of  simplicity,  the  ordinary  space 
and  imagine  to  have  established  a  system  of  general  curvilinear  coordi¬ 
nates  y1,  y2,  y3,  such  that  we  have  a  one-to-one  correspondence  between 
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A 

the  sets  y  *  (y1,  y2,  y3)  and  the  points  P  of  space.  It  is  well 

r 

known  that  the  equations  y  -  constant  give  us  the  three  families  of 

coordinate  surfaces^  each  surface  will  be  specified  by  one  and  only  one 

r 

value  of  the  parameter  y  .  The  intersection  of  the  surfaces  of  the 

three  families  give  us  the  "congruence"  of  lines,  the  coordinate 

/  r  r+lN 

lines,  each  of  which  is  specified  by  a  pair  of  coordinates  (y  ,y  ). 

t+2 

Along  each  coordinate  line  only  the  third  coordinate  y  may  vary. 

Let  us  consider,  furthermore,  a  point  P(y£,y*,yo)  *  P(yg),  ^ 
at  this  point  the  three  vectors 


T  -22. 
<  ’  1  < 


(III-28) 


These  vectors  have  obviously  the  directions  of  the  coordinate 

lines  through  the  point,  in  the  order.  Thus,  for  instance,  vx  is 

directed  along  the  intersection  of  the  surfaces  y2  «  y2  and  y3  m  y3  . 

The  direction  of  them  is  that  of  the  increasing  y’ s,  and  the  length 

ds 

is  given  by  the  ratio  — r  of  an  infinitesimal  displacement  along  the 

dy1 

line  and  the  increment  of  the  corresponding  parameter. 

Thus  if  a  general  system  of  curvilinear  coordinates  is  given  in 
space,  we  are  able  to  define  in  a  veiy  simple  and  natural  way  a  funda¬ 
mental  system  of  vectors,  one  for  each  point.  It  is  obvious  that  the 
systems  so  established  are  different  each  from  the  other  not  only  as 
their  orientation  is  concerned,  but  also  in  their  interior  structure. 

By  this  way,  we  are  now  able  to  define  in  each  point  the  components 
of  our  absolute  entities,  i.e.,  tensors,  by  simply  applying  the  theories 
stated  above  for  rectilinear  coordinates.  The  only  difference  is  that 
the  same  entity  (a  vector,  for  instance,)  will  not  have  the  same 
components  at  different  points  of  space,  as  was  the  case  in  rectilinear 


coordinates. 
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The  further  question  is  to  see  how  it  is  possible  to  connect 
the  different  fundamental  systems  of  vectors  we  have  established  at 
each  point,  each  to  another,  so  as  to  make  it  possible  to'  compare  the 
same  absolute  entity  at  different  points,  although  its  components  are 
different  from  point  to  point. 

Before  meeting  this  problem  we  may  observe  that  we  may  obviously 

write 

dP  «  ^dy*  4  y3  -  ^dy1  -  v±<fyj  ,  (111-29) 

and  thus  dy^  are  the  contra variant  components  of  the  elementary  dis¬ 
placement  dP.  (In  the  case  of  rectilinear  coordinates  yi  would  be 
the  contravariant  components  of  the  positional  vector  P-O;  this  is 
no  longer  true  in  general  coordinates,  and  is  replaced  by  the  foregoing 
differential  relation. ) 

Also  we  have 

dP.dP  -  ds2  »  v-^.VjC^dyj  »  g^dy^-dy^  ,  (III-30) 

having  written  for  v^.Vj  .  g-jj  is  therefore  our  fundamental 
metric  tensor.  It  defines  completely  the  internal  structure  of  the 
fundamental  system  of  vector*  (v^)  at  each  point.  As  we  can  see,  it 
depends  upon  the  position  of  the  point  considered  (in  the  case  of 
cartesian  rectilinear  coordinates  it  was  constant  for  the  whole  space). 

3.6  Connection  Coefficients 

Let  us  consider  a  point  PCy1")  of  space,  to  which  the  fundamental 

system  of  vectors  (vr)  -  ^r  13  related,  and  an  infinitely  near  point 

P'  (yr  4  dyr)  ■  P  4  dP  (dP  -  vj  dy^)  to  which  the  fundamental  system  of 

vectors  (v* )  -  — —  may  be  related.  *e  may  put 
r  Byr 

vi  -  Vi  4  dVi  (in-31) 

and  we  want  to  evaluate  d?i  . 
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As  a  first  step,  and  as  far  as  dv-j_  is  a  vector  attached  to  P, 
we  may  write 

dv±  »  Oh  v  ,  (III-32) 

x  i  k 

k 

( 0  .  being  thus  the  contravariant  components  of  dv^  at  P. 

Moreover,  u)  depends  upon  the  displacement  dp,  i.e.,  upon  its 
i 

(contravariant)  components  dy1-1 j  and  so  far  as  we  may  confine  ourselves 
to  a  first  approximation  (we  neglect  infinitesimals  of  higher  order) 

jc  1 

we  my  always  write  ^  as  a  linear  combination  of  dyJ  : 


k  rk  j 

«>i  *  /y  <iy  ,  (III-33) 


and  therefore 


dv  -  rk 

i  /  ij  k 


(111-34) 


r*k 

The  /  are  systems  of  the  third  order,  called  connection 

'  ij 


coefficients.  Their  knowledge  gives  us,  in  fact,  the  means  of  connect¬ 
ing  our  fundamental  system  of  vectors  to  the  neighboring. 

The  connection  coefficients  are  not  tensors,  as  may  be  easily 
understood  considering  that  a  tensor  is  essentially  an  entity  attached 
to  one  point  of  space  only,  whereas  the  p s  are  attached  to  the  point 
and  to  the  neighboring. 

We  may  find  an  algebraic  expression  for  the  connection  coefficients 
if  we  only  consider  that 


c)v. 


dv. 


3-V3  , 


(III -35) 


and  thus 


(111-36) 


if  w®  multiply  this  equation  by  v  w@  get 


Bv.  r  k  f>  ^ 


K  r  n  kK  r-K 

•7  *  H  \  “  /  •  •  » 

'  1 J  h  '  ij 


(III-37) 


which  gives  us  the  requested  formula. 

_  _k  *  k 

Let  us  remember  now  that  v  .v  ■  o  and  differentiate 

i  i 

partially  with  respect  to  y**: 


B  v.  _k  d  v  — 

_ i  -.v  +  *  v  *  0  j 

V  Vi' 

we  may  therefore  write  the  following  expressions  for  the  f~% a: 


rk  _k  V  _  3vi  _k  cfiF*  _ 

/y  "  V  -v  “  “  V  *  Vi  ”  V  *  v  “  ‘  VJ  “ 


•3  P  _k_  rK 
W‘v  /ji 


(111-38) 


The  p 3  are  also  called,  if  expressed  algebraically, 
Christoffel  symbols  of  the  second  kind.  The  Christoffel  symbols  of 
the  first  kind  are  expressed  by  the  following  relations 

,  .  r-k  dv±  _  Bv, 

UJ.h)  -  ghk/  -  (Ji, h)  "  3^-  'i  -  •’fc 


(III-39) 


Both  the  symbols  of  Christoffel  of  the  first  and  of  the  second 
kind  may  easily  be  expressed  in  terms  of  the  fundamental  tensor  g 

ij 

and  its  first  partial  derivatives.  But  we  are  not  interested  in 
developing  these  formula®,  that  are  to  be  found  in  each  book  dealing 
with  tensor  calculus. 


H  A  T)0  ^F,  c2— 149 


39 


We  may  only  give  a  formula  that  will  be  useful  in  the  following 
and  which  establishes  a  simple  relation  between  the  Christoffel  symbols 
of  the  second  kind  and  the  value  D  ■  Jg^jj  of  the  determinant  of  the 
fundamental  tensor. 

We  have 

dyk  ■  9yk  ■  Tj  *  Ti  -  •  (ik,j)  *  (*,i)  , 

and  moreover,  by  contracted  multiplication  by  g^: 


gij  Mil  =  r1  +  ra  »  2  r 

g  I  i*  I*  'rk  • 

3g  *  * 

We  now  observe  that  g1^  might  be  written  as  follows: 

ij  1  Gij 

g  3^*  IT  u  2>j*  ’ 

G  being,  by  definition  of  g^~\  the  algebraic  complement  of 

gii  in  lgiil  (Gld  m  D  glJ)* 


m  _  _ij  ^Si^ 

But  G  3yk  is  no thing  else  than  the  sum  of  three  determinants 
obtained  from  g-^ j  by  substituting  one  of  the  rows  with  the  derivatives 
of  the  corresponding  elements,  and  therefore,  as  it  is  well  known, 
this  development  gives  us  the  derivative  of  D: 

i  GiJ  !f±i  .  1  ®  .  31e  p 

^  Sjrk  D  3yjk  * 


And  finally  we  get  the  requested  formula: 


fr 


r  »  i-  9  lg  D  ,  3  lg  vt> 

rk  2  3yk  3yk 


(111-iiO) 
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In  a  three  dimensional  space  the  number  of  Christoffel  symbols 
of  each  kind  is  27  (  -  3*),  but  only  1.8  are  distinct,  owing  to  the  fact 
that  the  symbols  are  symmetrical  with  respect  to  the  two  lower  indices. 

We  will  still  observe  that  by  multiplying  dP  -  vr  dyr  by  vs  ,  we  get 

7*. dP  *  v  .  ^ dyr«  SSdyr  -  dy®  , 
r  r 

and  therefore,  by  the  definition  of  the  gradient,  we  may  write 

v®  *  grad  y8  ,  (111-1*1) 

i.e.,  the  reciprocal  vectors  of  each  fundamental  system  are  given  by 
the  gradients  of  the  parametric  coordinates  at  that  point. 


3.7  Absolute  Differentiation 

Let  us  have  a  field  of  vectors  u(P)  ■  u  vr  ,  each  vector  being 

r 

attached  to  the  corresponding  point  P;  w®  may  compare  Tt(P)  to  the 
neighboring  vector  u(P+dP)  ■  u(P)  +  du  .  As  a  particular  case»  we 
may  have  identically  dtt  -  0,  i.e.,  the  field  of  vectors  would  be 
formed  by  equipollent  vectors. 

From  u  *  u  v1*  we  have  by  differentiation  (interchanging  some 
r 

indices) 


cm 


dUj.  vr  + 


dvr 


(111-1*2) 


3U]r  _k  h 

and  not  du  *  v  dy  as  would  be  in  the  case  of  cartesian 
rectilinear  coordinates. 

r-r 

The  presence  of  the  term  is  due  to  the  f&ct  that  in 

in  passing  from  P  to  P  +  dP  not  only  u  has  varied,  but  also  the 
reference  system  of  vectors. 
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We  put 


Ifcl-  \  I  hk  “  Vh 


( 111-1*3) 


(read  u^  derived  h)  and  call  u^^  th®  ccvariant  derivative  of 
;  we  also  may  call 

)djrh 


(m-i*i*) 


the  absolute  differential  of  Uj.  . 

In  a  quite  similar  way  we  may  obtain  for  contravari&nt  components 
and  for  double  tensors  the  following  formulae 


k  W 

u  A  -  S3 


r  *~k 
U  Irh.  ’ 


*  “'A  *  *3^:  "  ahJ  rk[~  *jk  /"a  »  8tc* 


( 111-1*5) 


Going  back  to  our  first  differential  formula,  it  is  very  worthy 
to  be  noted  that  it  may  be  written  as  follows s 

d W  -  u^v  dy  ,  (111-1*6) 

and  that  this  formula  is  derived  from  the  expression  u  -  u^v  by 
the  ordinary  rule  of  differentiation  in  cartesian  rectilinear  coordi¬ 
nates,  only  substituting  the  covariant  derivative  u  ,  in  place  of 

8uk  k/h 

the  partial  derivative  :JyE  *  also  ®a/»  therefore*  that  the 

fundamental  vectors  behave  as  constants  with  respect  to  the  covariant 
(tensor)  differentiation. 

This  is  a  way  of  expressing  the  fundamental  theorem  of  Ricci. 

We  may  now  consider  the  particular  case  where  dTE  ■  0  identically, 
i.e.,  the  case  of  a  parallel  (equipollent)  field  of  vectors.  In  this 


case  we  have 


and  therefore 


u  ■  0  identically, 
k/h 


^k  r-r 

%S-Ur/hk 


(111-1*7) 
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3.  a 


If  w©  are  given  an  arbitrary  curve  (T)  in  parametric  fora 
y*  m  y“(s),  a  being  the  arc  sssasured  on  the  curve,  we  may  put  for  the 
sake  of  simplicity  »  X*.  Thus,  if  we  call  T  the  tangent  unit 


vector  to  our  curve  at  P,  we  have 

r  dP  ,i_ 

t  ■  T*  *  ~™  -  X  v.  | 
da  da  i 


(III-W) 


the  X  *»  are  therefore  the  contravariant  components  of  the  tangent 
unit  vector  to  a  curve. 

As  it  is  well  known,  Frenet's  formula  gives  us 
df  ui 


~  ■  Kn  *  KvSr  , 
d8  i 


(III-U9) 


K.  being  the  first  curvature  (or  flexion)  of  (r),  and 
K  »  y 'T_  the  unit  vector  of  its  principal  normal.  The  same  formula 
may  be  written  in  tensor  form  as  follows? 

XVr  Xl"  "  Kpl  <becau*e  ^  “  §r  jjj“)  .  (HI-50) 

If  K.  ■  0  identically,  our  curve  is  a  straight  line,  and  there- 


X*/r  **  -  0 

may  be  regarded  as  the  tensor  differential  equation  of  a  straight 
line  in  space  referred  to  general  curvilinear  coordinates.  We  may 
also  write  more  explicitly 

(  *  \k[~^  )\r  .  o  , 

'  9yr  '  i  r  ,  i 

and  moreover,  remembering  that  \T  *  ——  ■%■■■  , 

By.r  3yE*  ds  ds* 


£#  * 


(III-51) 


(III-52) 


This  is  the  final  form  of  the  differential  equation  of  a  straight 

line  which  we  will  use  in  the  following* 
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3»9  Some  Complements 


The  great  help  we  are  able  to  find  in  tensor  calculus  is 
essentially  due  to  the  fact  that  formulae  may  be  written  in  a  very 
compact  and  transparent  form  even  though  they  contain  a  great  number 
of  terms,  and  moreover  on  the  fact  that  the  following  rules  hold: 

(1)  A  zero  tensor  (the  components  of  which  are  zero)  in  a  given 
system  of  coordinates,  is  identically  zero  in  each  other 
system  of  coordinates  (in  fact,  the  transformation  between 
different  systems  is  linear  and  homogeneous), 

(2)  Two  tensors  (of  the  same  kind  and  order)  identical  in  one 
system  of  reference  are  identical  in  any  other  system. 

(3)  The  properties  of  symmetry  and  of  skew-symmetry  of  a  tensor 
are ir variant  for  ary  transformation  of  the  base  vectors  and 
indicate  therefore  an  intrinsic  property  of  the  entity 
represented. 

Therefore  a  tensor  equation  presents  an  absolute  character, 
whatever  may  be  the  system  of  coordinates  used  for  writing  it. 

We  will  at  once  show  some  consequences  of  the  above  rules. 

First  of  all,  we  may  very  easily  recognize  by  another  way  that 
the  Christoffel  symbols  are  not  tensors;  in  fact,  they  are  all  zero 
in  cartesian  coordinates,  and  this  is  no  longer  true  in  whatever 


other  system  of  coordinates. 


du 
dP  * 


Secondly,  let  us  consider  a  vector  u  *  u'Nr.  and  the  homography 
In  cartesian  coordinates  we  may  write,  as  we  have  already  seen. 


du 

dP 


dv.1 

du2 

du3 

ax1 

dx1 

du1 

3u* 

dx2 

cbc2 

dx2 

du2 

dx3 

dx 3 

dx3 

( III-53) 
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3.10  Curves  in  Space 

Let  us  have  a  curve  in  space  defined  by  its  parametric  equation 
P  «  P(s),  s  being  the  arc  of  the  curve.  We  immediately  have  that 


(III-57) 


is  the  tangent  unit  vector  to  the  curve. 

Considering  that  t.t  **  1,  we  have  t.—  *  0,  and  therefore 

—  t  ds 

t  and  ^  £  are  orthogonal.  We  put 


dt  d*P  n 
ds  dsz  P 


(II 1-58) 


TI  being  the  unit  vector  of  the  principal  normal  to  the  given  curve, 
and  p  being  the  radius  of  first  curvature  (or  flexion)  of  the  curve. 
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Let  us  next  call  b  the  unit  vector  normal  to  both  n  and 
t,  and  such  that  the  system  (t, n,b)  may  be  positive.  We  have  from 
b.b  -  1  and  fc. t  *  O: 

b  -  0  .  —  •  t  +  «  0,  and  therefore  also  ~  .  t  «  0. 

ds  ds  y  ds 

J'U  _  ___ 

That  means  that  —  is  perpendicular  to  both  b  and  t  and  is 
ds 

directed  therefore  along  n.  We  put 

db  _  n 
ds 

and  call  x  the  radius  of  second  curvature  (or  torsion)  of  the  curve. 
From  n.n  «  1,  n.F  »  0,  n.F  *  0  we  get  after  this  by 
differentiation,  taking  into  account  the  foregoing  formulae, 

~  .  n  -  0,  .  t  «  -  -L  —  .  F*  2- 

ds  ds  p  *  da  X  » 


and  therefore 


dn  t  b_ 
ds  "  ”  p  ~  T  * 


The  set  of  formulae  which  we  summarize 

d  t  _  n  dn  _  _  t_  _  b_  db  _  n 
ds  f>  ’  ds  ~  p  ~  Z  ’  ds  = 

are  known  as  Frenet* s  or  Serret1 s  formulae. 

Tne  plane  of  t  and  n  is  called  tangent  plane  to  the  curve 
the  plane  of  t  and  b  is  called  the  rectifying  plane;  and  the  plane 
of  n  and  b  is  called  the  normal  plane. 


(III-59) 


3.11  Curves  on  a  Surface 

Let  us  now  consider  a  curve  (  T  )  on  a  surface  Z;  we  will 
indicate  by  N  the  normal  unit  vector  to  Tt  and  by  n  the  principal 
normal  to  (T).  t  being  always  the  tangent  unit  vector  to  (  71) 
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(perpendicular  to  both  n  and  N),  we  may  also  consider  the  positive 
trihedron  (surface  trihedron)  (T,!T,E)  of  normal  unit  vectors,  and  put 

H.n  -  cos  ©  ,  B.b  *  cos  0  , 

n.3  ■  sin  6  ,  N.b  *  •  sin  0  , 

0  being  the  angle  of  obliquity  of  (  l1)  on  Z. 

Therefore; 


n  ■  N  cos  0  +  B  sin  0  , 

b  »  -N  sin  0  +  B  cos  0  , 

N  «  n  cos  6  -  b  sin  0  , 

B  ■  n  sin  0  +  E  cos  ©  . 


(III-60) 


By  substituting  these  fomulae  in  Frenet' s  formulae  we  easily 


dl  TT  .  1 


ds  fn  pg 


ds  Tg  fa  ’ 


(III-61) 


dB  .  J_  _  X 

da  Pg 

which  are  called  the  generalized  Frenet' s  formulae  for  the  surface 
trihedron;  and  where  we  have  put 


1  cos  Q 

Pnm  P 


(normal  curvature)  , 


sin  0 

P 


(tangential  or  geodesic  curvature)  ,  (III-62) 


—  +  —  (geodesic  torsion)  . 
x  ds 


For  geodesic  curves  on  the  surface  we  have  identically  n  ■  N, 


and  therefore  0  -  —  ■  0;  thus  for  geodesic  curves: 

ds 

1_  1_  X  n  x  -i- 

P  "  fa  9  Pg  U  *  Tg  "  Z  • 


( III-63) 
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3.12  Surfaces 

Let  us  consider  a  surface  Z  and  its  normal  N,  and  take  an  arbitrary 
displacement  dP  on  2.  If  we  give  the  point  P  the  displacement  dP, 

TT  would  become  N  +  dN;  and  from  N«N  ■  1  we  have  that  N.dN  *  0,  i.e., 
dS  is  a  tangential  vector. 


The  homography  (of  Burali-Forti) 

a.  € 
dP 


(111-61*) 


such  that  crdP  «  dN  gives  us  the  law  of  variation  of  N  in  any  direction 
and  defines  therefore  the  shape  of  Z  in  space. 

If  we  consider  the  parametric  form  of  2, 

p  -  P(y\  y2), 

y1  and  y2  being  two  parameters,  we  have 


dP  -  ~  dyl  +  ^  dy2  *  vi  dyl  +  v2  *  7a  dy°  "  (III-65) 

and  therefore 

2  2  a 

ds2  -  dP.dP  -  Edy1  +  2Fdy1  dy2  +  Gdy2  -  dya  dyP  , 

(first  fundamental  form  of  Z) 
having  put 

gja  -  E  -  va  .  va  ,  gX3  -  g2X  -  F  -  va  .  v2  -  v2  .  vx,  g23  -  G  -  V2 .  v2.( III-66) 

The  g^(a,P  =  1,2)  are  the  components  of  the  metric  (first) fundamental 
tensor  of  2  in  our  coordinates  ya. 

We  may  observe  that  F  -  \/EG  cos  0,  where  6  is  the  angle  of 
the  considered  coordinates  lines  on  2;  and  therefore 


cos  e  .  -E.  ,  sin  e  -  /— 2s .  ( 111-67} 

Considering  now  c%  we  may  put 


°  Td 


( HI-68) 
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and  call  b ap  the  second  fundamental  tensor  of  Z.  We  have 


rf«  ■  T3  *  -  !  <*  '  -  bac 


(III-69) 


and  moreover 


dN.dP  »  -  bap  dytt  dyP  (second  fundamental  form  of  Z). 
Finally  we  have 


(HI-70) 


dfi.dS  *  bfl^  v°  dyP.  vr  dyJ  *  gaT  ba  dyf  dy*  -  dyP  dy*  (III-71) 


(third  fundamental  form  of  Z) 


and  *  ga  3  b^  b^  is  called  the  third  fundamental  tensor  of  Z. 

We  may  observe  that  -dN.dP  -  N»d*P  because  d(N.dP)  - 
= dN.dP  +  TI.d*P  *  0  ,  and  also  that  from 

M.dP  -  N.7a  dya 

we  get  by  partial  differentiation  and  interchanging  the  indices  ct 


and  |3  s 


(HI-72) 


3ya  P  ^  3yaayP 


dyP  -  0  , 


and  therefore 


3N  _ 

£T'Va  =  -  b 


_  3N_  _ 

ajl  "  -yfi  ’  ~  ba£ 


( III-73) 


The  second  fundamental  tensor  is  therefore  symmetric  and  is 
usually  indicated  as  ba^  *  (D,D' ,D"  ).  Also,  it  follows  that  the 
homography  <T  is  a  dilatation,  owing  to  the  fact  tnat  the  property  of 
symmetry  of  tensors  is  invariant  with  respect  to  any  transformation. 
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We  remember  now  that 


2L.*  r' 

dya  aP 

and  differentiate  this  formula  covariantly  with  respect  to  y* 
Remembering  also  that  the  fundamental  base  vectors  (and  their 
reciprocals)  may  be  regarded  as  constants  with  respect  to  covariant 
differentiation,  we  get 


3N  d2  N  3N  r &  -p 

dya/p  ”  3ya3y^  3yS  ~  */>//* 


Interchanging  a  and  j3  ,  the  central  part  of  the  equation 
written  remains  unaltered  (on  account  of  the  symmetry  of  the  second 
partial  derivatives,  and  of  the  Christoffel  symbols  of  the  second  kind 
with  respect  to  the  lower  indices);  and  therefore 


bcc/>/p  *  hafi/p  * 


(III-71|) 


This  expresses  in  tensor  form  the  Mainardi-Codazzi  equations, 
which  give  the  integrability  conditions  which  the  second  fundamental 
tensor  of  a  surface  must  satisfy. 

Moreover,  if  we  take  a  displacement  dP  in  the  direction  of  the 
unit  vector  t  on  Z,  from 


we  immediately  have 

1 


(T  1. 1  *  - 


rt.B  - - 


<r  dF  _  ,  dh  „  _  _3_  _  _t_ 

ds  Tg  Pn  * 

that 

(normal  curvature  of  Z  in  the  direction  t  ) 
(geodesic  torsion  of  Z.  in  the  direction  T). 


(ixi-75) 


(1 11-76) 


The  united  directions  c%  and  c2  of  Zf  such  that 


h  9 


( III— 77) 


(which  always  exist  since  <f  is  a  dilatation)  give  us  the  principal 

directions  of  (T }  p1  and  p2  are  the  principal  radii  of  curvature 

of  Z,  Along  these  directions  the  geodesic  torsion  is  zero. 

Also  I,<r  gives  us  the  mean  curvature  “■  +  ~~  ,  and  I_<T 
1  Pi  2 


the  total  curvature 


of  Z. 


V.'AIX 


Th  52-1^9 


51 


SECTION  IV 


INTRINSIC  GEODESY 
(See  reference  i*i*,  51,  52,  53) 

k*1  Intrinsic  Geographic  Coordinates  on  a  Surface  (See  reference  26.  2?) 

If  we  take  a  constant  anticlockwise  (sinistrorsum)  trihedron  of 
orthogonal  ypit  vectors  (a^  ag,  *3)  and  a  variable  clockwise  (dextrorsum) 
trihedron  of  orthogonal  unit  vectors  (ix,  i2,  i3),  we  may  always  write 

11  *  “  sin  ♦  cos  X  .  ix  -  sin  $  sin  X.  a2  +  cos  ♦  .  a3  , 

12  *  “  sin  X.  ax  +  cos  X.  a2  , 

13  -  cos  *  cos  X  .  5^  ♦  cos  $  sin  X .  a3  +  sin  * .  a3  , 

(IV-D 

*o.  “  “  sin  ^  cos  X .  -  sin  X  .  i2  ♦  cos  #  cos  X .  i3  , 

I2  *  -  sin  sin  X  .  Ix  +  cos  X  .  I2  +  cos  $  sin  X .  i  , 

a3  "  cos  *  •  Ii  +  sin  <t .  I  , 

where  $  and  X  are  the  latitude  and  longitude  (respectively)  of  i3  with 
reference  to  a3  as  an  axis  (world's  axis  directed  towards  the  North) 
and  the  plane  (a2  t  a3)  as  fundamental  meridian^  and  we  understand 
that  ix  must  always  be  in  the  plane  (i3  ,  a3),  i.e.,  directed  towards 
the  North. 


By  differentiation  we  get  immediately 


di2  *  i3  sin  $.  dX 
di3  -  I1.  d» 


-  i2  sin  $.  dX  -  i  ,  d*  , 

-  I3  .  cos  #  .  dX  , 


( IV- 2) 


♦  i2  cos  *.  dX  , 


aind  therefore 


dl3  ,  di3  *  (di3)2  »  d*2  +  cos2$dX2  . 


(IV-3) 
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Let  us  now  take  a  surface  Z,  and  its  normal  I  at  a  point  P 
and  establish  a  correspondence  (Gauss's  spherical  correspondence)  such 
that  N  ■  i3  .  Under  a  certain  hypothesis  of  convexity  for  the  surface 
which  we  will  assume,  the  correspondence  will  be  one-to-one;  and  we  also 
may  therefore  attribute  to  P  the  coordinates  ♦  amd  X  of  I3.  le  will 
call  them  the  intrinsic  geographical  coordinates  of  P  on  Z.  In  the 
case  in  which  Z  is  the  geoid,  the  intrinsic  geographical  coordinates 
coincide  with  the  astronomical  coordinates. 

We  then  have 


fi’dP  -  dN  *  ij  d#  ♦  i2  cos  #  dX  ; 
ix  and  i3  give  us  the  North  and  East  directions  on  Z. 

If  we  now  consider  the  fundamental  vectors  on  2.  connected 


3P 

9X 


with  the  intrinsic  geographic  coordinates,  vx  «  ,  va 

(the  directions  along  which  the  longitude,  and  latitude 
do  not  vary,  i.e.  the  geodetic  meridians  and  parallels  on  Z),  we 
immediately  have  from  the  foregoing  formula®  for  our  second  fundamental 
tensor: 

,  n  3N  -  T  - 
‘’n  '  D  ‘  ""5* ‘ Ti  ‘  1»-Ti 

*>1.  "  D'  ■  -41-  •  -  ii-  ,  tfc. 


.  aff  _  _  _ 

,  b21  -  D  -  -  —  "  i2  *  cos  *  , 


D"  -  -  .  v„ 


ax 

3N 

ax 


_  _  (i V-ii) 

-  .  Ta  cos  *  , 


and  therefore 


ix  D  -  i2 


T  D' 


cos  f 


-  ^  »  -  ia  — 


(iv-5) 


C03  $ 
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If  we  multiply  these  relations  and  remember  that  gu 
immediately  get 


E  -  ya  . 


gxx  -  E  •  D* 


D»‘ 


cos3#  * 


D* 


■  «ax  ■  »  ■  »(»  *  35  )  , 


(IV-6) 


gaa  -  G  -  D»3 


cos8# 


Although  in  general  the  first  and  the  second  fundamental  forma 
are  not  related  by  equations  in  finite  terms,  in  this  case  (of  geographic 
intrinsic  coordinates)  they  are,  owing  to  the  particular  spatial  meaning 
of  the  coordinates  themselves. 

As  may  be  seen,  the  total  curvature  of  £  is  given  by 
„  1  cos8# 


Pip* 


DDm  -  D‘ 


( IV— 7) 


Two  vectors  a  and  T>  such  that 

a  .  <rb  -  0  ( IV- 8) 

and  said  to  be  conjugate  with  respect  to  cr  f  it  follows  from  the  fact 
that  cr  »  Ker  (because  <r  is  a  dilatation),  that 

a  .  <r  b  -  b .  era  , 

and  therefore  the  property  of  being  conjugate  is  reciprocal. 

We  have  from  the  foregoing  that 
dS  ■  ix  d#  +  i2  cos  *dX  *  <rdP  -  «r  vx  d#  +  cr  va  dX  , 
and  therefore 


o-vx  .i3  -  ^  .ia  -  vx  ,o-ra  -  0  , 

<rva  .±x  -  I2  .1  cos  #  ■  v  .  cri,  -  0  , 


(IV-9) 
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i.e.,  the  directions  of  vx  end  v3  are  conjugate  respectively  to  the 
directions  of  i2  and  ix.  This  leads  to  the  theorem  first  formulated 
by  Pizzetti:  the  direction  of  the  meridian  (parallel)  on  ,a  given 
surface  is  conjugate  to  the  East  (North)  direction  with  respect  to  the 
homography  of  Burali-Forti  (or  also  with  respect  to  the  indicatrix  of 
Dupin).  (See  reference  31r  32) 

We  also  may  observe  that 

—  -  / —  P  D» 

v,  ,  i,  -  vE  cos  p.  »  -D  .  cos  u  »  -  ,  tg  p  -  ~ - 

vE  D  cos  * 


—  —  r-  V”  D« 

v2  .  i_  -  vQ  sin  *  - - —  ,  tg-rr  - - 

3  ^  cos  #  D*  cos  ♦ 


(IV-10) 


give  us  the  aximuths  jjl  and  ir  of  the  meridian  and  of  the  parallel, 
respectively,  and  that 

r—  D«  D‘  D" 

?X  ,  v2  -  cos  9  -  D'(D  +  - )  ,  cos  8  ■  -=  (D  ♦  -  )  (IV-ll) 

cos  v  VEG  cos  ♦ 

give  us  the  angle  0  between  the  meridian  and  the  parallel. 
ii.2  Curves  on  a  Surface  Referred  to  Geographic  Coordinates 


If  we  remember  the  generalized  formulae  of  Frenet  and  put 
t  ■  ±x  cos  A  +  ia  sin  A  (A  •  azimuth  of  T) 


—  ha  —  HA  dli  dl_ 

-  -  i*  sin  A  —  ♦  i  Cos  A  ~  ♦  -r-~  cos  A  ♦  — *  sin  A  , 
ds  *  ds  ds  da 


(IV-12) 


we  immediately  get  the  important  formulae  valid  for  an  arbitrary  curve 


on  Z  : 


1  .  _  dX  .  dA 

— -  »  -  sin  ♦  —  ♦  —  . 

Pg  ds  ds 


1  *  ,  *  dX  » 

-  —  “  cos  +  sin  A  —  ♦  cos  A  — 

P-n  ds  ds 


( IV-13) 
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For  a  geodesic  we  have  therefore 


ln  dA  .  ,  dX 

fig  da  ds 


(IV-lJ*) 


Moreover,  we  may  write 

-  , d—  .a  dP  m  d_P_  f 

1  “  X  Ya  •  X  3ytt  3y«  ds  *  (3 

and  therefore  X°  »  are  the  contravariant  components  of  the 

ds 

tangent  unit  vector  t  to  a  curve  P  *  P(s)  on  Z  .  If  this  curve 
would  be  a  geodesic,  the  ya,s  would  satisfy  to  an  equation  similar  to 
that  we  found  in  3*8,  i.e., 


(iv-15) 


fr- 


(IV-16) 


This  formula  may  be  differentiated  successively;  remembering 


that,  for  instance. 


s&Lsi! 

a*  dsa 


we  easily  get 


d.»  (2  V  £»  ‘a/  )xxx  * 


(IV- 1?) 


These  developments  allow  us  to  generalize  for  a  geodesic  line 
drawn  on  a  surface,  Legendre’s  series  for  the  computation  of  geographic 
coordinates  from  polar  geodesic  coordinates,  (see  reference  hi)  i;e.. 


ya  -  yj  +  * 


X«*£  (*£)  +»-  (£& 

0  2!  \  ds  Jo  31  \  ds 


+  •  *  *.  (IV-18) 


For  the  azimuth  we  get  in  a  similar  way  the  following  expansion: 


'•-‘MtVS®).---. 


(IV-19) 
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The  c ontravar iant  c onponents  k  are  easily  expressed  in  terms 
of  the  aximuth  A  as  follows: 


V1  i  t 


•  1, 


coe  A  «  A3^  •  1^+  Aava»  ix 
-  DX1  -  D* X2  , 


t  •  ia  *  sin  A  *  - 


D*  X1  ♦  D"Aa 


and  therefore 


(-  D*1  cos  A  +  D*  cos  #  sin  A)  , 
(B*  cos  A  -  D  cos  $  sin  A)  , 


( IV-20) 


A  -  DS"  -  D*2  . 

Formulae  perfectly  similar  to  these  may  also  be  established 
for  a  curve  on2  ,  provided  the  geodesic  curvature  be  known  as  a 
function  of  the  arc  s. 


U.3  Intrinsic  Geodesy  in  Space;  Absolute  Formulae.  (See  references 

16,  1^) 

Let  us  consider  now  the  vector  g  of  gravity  in  space  and  put 

I  -  -  gn  *  (IV-21) 

as  it  is  very  well  known,  n  is  the  normal  unit  vector  to  each  equipo- 
tential  surface  of  gravity,  and  g  is  the  intensity  of  gravity.  We 
take  g  and  fi  in  opposite  directions. 

If  we  differentiate  the  foregoing  formula  with  respect  to  P  , 
we  obtain  a  homography  (Ebtv&s's  homography) 

^  "  dl"  &  grSd  ¥*  ( IV-22) 


WADC  TR  5 2-149 


57 


(IV- 23) 


f  b#ing  obviously  tbs  potential  of  gravity^  w®  also  will  writ® 


and  call  cr  the  generalized  holography  of  Burali-Forti.  (generalised* 


because 


cr  is  no  longer  related  to  a  single  surface,  but  to 


family  of  surfaces,  the  equipotential  surfaces). 

We  immediately  get 

dg  dn  -  dg  _  -d£ 

dP  dP  dP  6  dP 


(IV-2U) 


where  the  meaning  of  —  is  given  by  the  formula: 

dP 

dP  -  dg  . 
dP 

The  operator  ^  is  therefore  identical  to  the  operator  grad  g, 
dP 

an  hyperhomograply  operating  between  vectors  and  scalars. 


dP 


grad  g  ; 


grad  g  is  called  the  gradient  of  the  intensity  of  gravity 
(gradient  of  gravity). 


Our  formula  may  now  be  written,  remembering  the  dyads,  as 


follows: 


and  therefore  also 


ZJ  -  -*  g  cr  -  H(grad  g,  n)  , 


(IV- 25) 


<U1UL  UJWA  OJLUiO 

*  H(grad  g,  U)J  . 

The  relations  between  ZJ  and  cr  now  established  give  us  in 
the  simplest  way  the  relation  between  the  mechanical  ( ZJ )  and  the 
geometrical  (  cr  )  structure  of  the  gravity  field  if  we  confine  our 
attention  to  its  characteristics  of  the  second  order. 


( 17-26) 


(IV- 27) 
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( 37-28) 


Ths  g ra  vime  trie  gradient  ssa^r  ia^sdlately  b#  @xpr©s*®d  by 
aeans  of  T<f§  we  observe,  therefor®,  that  for  asy  taro  vectors  x  sjsi 
f  we  hftv®9  r6»s^b«?iEg  Jacobi*  a  theorem  ( 11-31)  and  ( 11—35)  , 

grad(x.y)  .dP  -  d(x.y)  -  f.dx  +  x.dy  *  F  •  ^  dP  ♦  x  .  ^  “ 

®  K  *  d?  *  K^x.dP  -  (K  ?f  ♦  K  ^  X)  •  dP 
dP  dP  ar  ay 

(K  —  is  the  conjugate  of  ~  ).  If  *  and  y  are  the  gradients  of 
dP  dr 

two  scalar  functions,  then 


dP  dP 


and  we  have  in  this  case 


dP  dP 

(17-29) 

Sy  .fix  . 
dP  dP 

(17-30) 

Thus 


grad(g.g)  -  grad  g*  -  2g  grad  g  ■  2  ^  f  -  2Ug  , 


and  therefore 


grad  g  ■  —  ZJ  y  •  —  ZlZ  n  • 

g 


We  also  have  in  consequence 


-i.ju'-  H(irs,a)J  . 


We  may  observe  that  the  homography <T  is  a  dilatation  if 
considered  on  a  surface  (as  already  seen),  but  it  is  no  longer  a 
dilatation  in  space.  We  have  in  fact 

%  <r  -  -  ” H(H,  grad  g)J  /  <*"  , 

because  K  U  -  V' j  and,  therefore, 

K  <T  h  .  0  , 

and  n  is  the  null  direction  of  Kf  (but  not  of  <r  )  % 


(17-31) 


(37-32) 


(17-33) 

(I7-3U) 
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The  vector  grad  g  is  always  perpendicular  to  the  surfaces 
g  ■  constant*  which  are  called  isogravitational  surfaces. 

Since  g  is  a  gradient,  we  have  furthermore 

rot  g  -  0  -  -  rot  g  n  *  -  g  rot  n  +  n  *  grad  g; 

(IV-35) 

rot  n  -  —  n  x.  grad  g. 

E 

We  may,  moreover,  observe  that 

rn  *  -  —  \xjn  +  H(grad  g,n)j 

g  1  J  (17-36) 

grad  g  -  (grad  g  .  n)  nj-  ~  grad£  g  , 

gradE  g  being  the  component  vector  of  grad  g  tangent  to  the 
equipotential  surface  £  (also  called  the  surface  gradient  of  g ,  or  the 
horizontal  gradient  of  gravity). 

Owing  to  the  fact  that  H  x  n  *  0,  we  may  therefore  also  write 


rot  n  »  —  n  x  grad  r  g  »  Kx  crn  . 

g  r 

The  vector  crfi  *  —71  gives  us  obviously  the  derivative  of 
dP 

the  normal  unit  vector  H  along  the  line  of  force;  its  modulus  is 

therefore  -  F  (curvature  of  the  line  of  force)  and  its  direction 

that  of  the  principal  normal  N'  to  the  line  of  force  itself.  If  we 

call  furthermore  5*  the  binormal  to  the  line  of  force,  we  will  have 

-  _  B' 

n  x  <r  n  ■  rot  n  ■  — 

R  * 

7  -  7  &rad  I  6  •  »'  • 

Considering  now  the  divergence  of  n  ,  we  would  have 

diw  n  -  -  div  —  div  g  -  g  .  grad  — 

g  g  g 

.  (l*7rf«r  -  2o>2  -  grad  g  .  n  ) 
g 

(f  *  constant  of  gravitation,  &  -  density  of  matter). 


(17-37) 


(IV-38) 


(17 -39) 
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The  divergence  of  n  has  a  precise  meaning;  let  us  consider 
the  homography 

crz  -  cr-  H(n,  ffn)  -<r[l  -  H(n,n)]  ;  (IV-1*0) 

we  have 

n  *  0  }  a^t  ■  crt  (t  ■  tangential  vector);  {TJ-kl) 

i.  e.,  crt  operates  on  tangential  vectors  like  cr  and  gives  a  null  vector 
if  operating  on  nonnal  vectors  to  £  .  We  call  therefore  crT  a  surface 
holography;  w©  immediately  see  that  it  coincides  with  the  homography 
of  Burali-Forti  for  £  which  vre  had  already  considered. 

He  also  have 

Ix<Tt  -  Ix«*“  Jx  H(n,  cr  n)  -  If  -  n  .  tr  n  -  l^er  ( IV-U2) 

since  n  and  cr  n  are  orthogonal  vectors.  Ix cr^  equals  now  the  mean 
11  c- 

curvature  H  *  —  +  —  of  Z.  therefore 

>°a 

H  -  div  IT  -  —  (i*irf$  -  2o>*  -  grad  g  .  E)  (VI-k3) 

& 

This  is  the  celebrated  formula  of  Bruns,  which  may  also  be 
written  as  follows 

grad  g  .  h  «  H  -  “-(2  po2  -  U  ir  f  S' )  .  ( IV— UU) 

o 

The  first  member  is  often  called  (improperly)  the  vertical 
gradient  of  gravity. 

U.1+  Local  Astronomical  Reference  Systems 

Let  us  consider  at  a  point  P  the  local  system  of  orthogonal 
unit  vectors  (i1>i2,i3),  at  which  as  usual  i3  is  directed  (like  n) 
towards  the  zenith  of  the  equipotential  surface  2  at  P,  i3  towards 
the  North,  and  i2  towards  the  East.  Along  the  axes  (1^,  i2,  i3)  we 
imagine  to  measure  the  cartesian  isometrical  coordinates  x1,  x2,  x3. 
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We  will  furthermore  put 

3*  m  w  a2  w 

3xr  *  r  *  dxrdx8 


3Wr  _  £5 

3  x8  3x 


8.  ■  IT  =  w 

r  rs  sr  » 


and  also  we  will  write 


( IV-^5) 


(IV-1*6) 


We  immediately  have  that  g  *  -  W3  ,  and 
g  «  grad  Kf  *  -  g  i3  , 

II  Wax  Wl*  W13 


XJ-  m  .  §&  «  w  w  w 

dP  dP  21  22  2311  * 


WWW 

"31  "32  "33 


( IV— 1*7) 


wrs  *  V  ±T  •  ig  *  ^  ^  •  *r  » 

tJ"!  -  w  I  +  f  7  ♦  W  I  , 

1  H  1  12  2  13  3  > 

Wla  ■  »aA  *  *aaTa  *  *, . 

'  WaA  *  Wa»ra  *  "a/a  • 

The  formulae  written  justify  the  name  of  Eotvos  homography 
given  toU  ,  as  Eotvos  was  the  first  geodesist  who  considered  sys tematically 
the  second  derivatives  of  the  potential. 

Furthermore 

div  g  -  WX1  +  W22  +  W33  -  2 oj2  -  liW  f5, 
grad  g  =  -Ui3  -  -  w3li;i  -  "  w33i3  » 


r  w 

11  12 


21  22 


( IV-48) 


r  w  0 

31  32 


and  therefore  we  get 

H  -  I xtr  m  Wli  +  w22  s  k 
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wnw2a  ~  A.z 


(I 7-h9) 


k,5  General  Intrinsic  Coordinates 


The  directions  forming  the  local  astronomical  trihedron 
(ij_,  i2,  i3)  have  already  an  intrinsic  character  and  may  therefore  be 
materialized  by  actual  physical  measurements  proper  to  Astronomy  and 
Geodesy;  but  this  does  not  imply  that  the  same  trihedrons  are  the 
principal  trihedrons  of  a  triple  congruence  of  coordinate  lines.  In 
fact,  we  would  not  be  able  to  find  three  congruences  of  curves  inter¬ 
secting  in  a  point  and  having  at  this  point  the  given  directions  of 
(i-U  ^-2>  ^-3)  • 

We  will  overlook  the  proof  of  this  general  proposition  and 
will  confine  ourselves  to  observing  that  if  we  consider  any  equipotential 
surface  and  draw  on  it  in  a  point  the  North  and  East  directions,  the 
normal  i3  along  these  directions  would  not  cut  in  general  the  infinitely 
near  surface  along  North  and  East  directions. 

The  astronomical  local  systems  of  vectors  are  therefore  not 
suitable  for  defining  a  general  system  of  coordinates  in  space. 

We  will  instead  chose  as  general  coordinates  in  space  the 
following : 

(1)  the  intrinsic  (astronomical)  latitude  on  each  equipotential 
surface,  *  or  y1  ; 

(2)  the  intrinsic  (astronomical)  longitude  on  each  equipotential 
surface,  A  or  y2  ; 

(3)  the  potential  W  or  y3  . 

The  three  families  of  coordinate  surfaces  are  therefore 

(1)  the  surfaces  of  equal  latitude  , 

(2)  the  surfaces  of  equal  longitude  , 

(3)  the  equipotential  surfaces. 
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We  may  observe  that  the  parameter  I  (the  potential)  is 
measured,  with  the  exception  of  an  inessential  additive  quantity  and 
a  factor,  by  the  dynamic  height. 

The  three  congruences  of  coordinate  curves  are 

1.  the  geoidal  meridians  (along  which  only  ♦  is  varying)  , 

2.  the  geoidal  parallels  (along  which  only  X  is  varying)  , 

3.  the  isozenithal  (or  isovertical)  lines  (along  wnich  only 
W  is  varying)  . 


The  three  fundamental  vectors  are  therefore 

ap  gp  _  ap  a  P  -  3P  gj* 

^  “  at  ‘  a7  *  v*  ’  ax  "  iy*  ’  ^  "  ***  sy3 


(iv-5o) 


The  direction  of  v  is  the  null  direction  of  &  ;  in  fact, 
<fdP  *  dn,  and  therefore,  since  v3  is  the  isozenithal  direction, 

v3  -  0  . 

If  we  would  attempt  to  express  the  fundamental  vectors  (vr  ) 
in  terms  of  the  local  astronomical  unit  vectors  (Ir),  we  only  want  to 
remember  that  v^  m  <rv2  ■  i3  cos  ♦  ,  tf"v3  »  0,  and  moreover,  that 

3P  -  3P  T  =r  -  1 

grad  —  -  g.  •  ^  -  g  i3  •  v3 

because 

grad  .  dP  »  dW  . 
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and  remember  that 


^  ix  -  ’ll1!  +  "x^a  » 


^ia  “  "ax1!  +  waa1a  » 


^ia  “  +  W3213  » 


we  will  have 

s-v,  =  Ix  *  (anwu  +  ai2W2i^i  +  (axlV**  *  *“W“)Ia 
crV2  *  i2cos  $  -  (a3lwll  +  a22w2l)i1  +  (aglwl2  +  a22w22)i2  , 

O-V  -  0  -  (a3iwil  +  a32W21  +  a33W3l)i:i  +  (a3lWia  +  a32W22  +  a33W32)j 


and,  therefore,  immediately 

w22 

axx  K 


D  }  aX2 


wi2  -5* 


cos  4? 


wl2 


COS  ♦ 


wlX  cos  ♦  -  D" 


21 


K 


-  D*  ;  a22 


cos  4f 


cos2  * 


K  *  wllw22  -  ^12  "  ^  i  A  *  D  D"  -  D1  , 


and,  moreover. 


a3lWH 

+ 

a32W21 

+ 

a33W31 

a31Wl  2 

+ 

a32W22 

+ 

*33*38 

- 

a33e 

and  thus 


W13W22  “  W23W12 

gD 


W23W11  ”  W13W12 


32 


*33 


gD 


-Hi  > 


-  h2  , 


1_ 

g 
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(IV-52) 


( iv-53) 


(iv-55) 


(IV-56) 


We  thus  have  finally: 


v  ■  -  D  I  -  T  , 

i  1  cosi  1  ’ 


-  D'ix  * 


D" 


cos  ♦ 


l2  » 


(IV-57) 


v  *  ~  H,  i.  -  HJL  -  —  i 


11  2  2  g  3 


The  meaning  of  and  Hg  may  be  seen  by  observing  that 
o- 7  -  0  -  -  -  K2<rla  -  J°-i3  , 


and  that 


_  1 

cr  i3  -  <r  n  *  —  grad£  g 


Therefore 


ix  ♦  Ha<rla  *  8rads  g 

_  1  _ 
and  multiplying  by  vx  and  v2  ,  remembering  that  grad£  —  •  v±  *  g^f 

and,  furthermore,  that  <r  i^v.^  **  .0vx  *  1,  <rix .  va  “crig.  vx  *  0, 

<rva .  ia  ■  cos  *  , 


we  get 


«i 


d  4> 


H. 


*  2  cos  i  H  * 


(IV-58) 


and,  therefore. 


vi  *  "  D  ^  " 


a?  T 


cos  $  2  * 


D*i - — —  i  , 

1  cos  i  * 


(IV-59) 


•^1  -V  1 

fri  _ _ i_ .ii  , 

£  ♦  1  cos  $  3  X  2  g  5 
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and  also 


5is  *  21  c 

4  ^  ^ 


—  D'  cos  $  —  D  cos  4f 


Ti  ~ 


» 


Sh  cos  $  —  g  k  cos  4>  — 

■■■'  *  . — "  v  +  v  — ' —  v 

A  1  A  2 


-  g  v 


3  * 


having  put  for  simplicity 

h  «  — i-  _  D.H  .  -  -  —1—  ( D"  i-  -  D*  Ar  ~)  > 

cos  ♦  2  g  cos  ♦  cos  ♦  3$  g  g 


*23  A 


k  -  DH  - 

3  cos  $  g  cos  4>  cos 


1  /n  3  1  3  1  \ 

j  (  3X  g  g^g) 


J4..6  Reciprocal  Vac  tors;  Components 

We  will  now  consider  the  reciprocal  vectors  v1*  t  given  by 

relations  vr*vs  *  S  ;  we  get  by  easy  computation 
8  s 

D"  -  .  D'  cos  *  t  g  h  cos  ♦  — 

v1  -  grad  t>  -  -  —  ix  ♦  -  i*  ♦  - * -  i3  » 

_o  .  .  n«  -  D  cos  §  -  g  k  cos  $  T 

7*  .  grad  X  -  2-  ^  -  — —  t  .  1 


v3  -  grad  W  - 


and  also 


-  D  v1 


v1 


2  cos  # 


al 

-  D'v2 - §-  v3 

a  *  » 


D" 


cos  # 


v2  - 


iLt* 

cos  #  3X  * 


.if 

g 


(IV-60) 


(IV-61) 


the 


(IV-62) 


(IV-63) 
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The  contravariant  components  dyr  *  (d4,  dX,  dW)  of  a  dis¬ 
placement  dP  -  i^dx1  +  i2dx2  +  i3dx3  are  related  to  the  local 
cartesian  coordinates  dxr  by  the  formulae 

d*  .  dy3  .  -  E.  dx3  .  *  6  h  dx3  , 

A  A  A 

,  D  C0S  *  ,*„2  ^  g  k  COS  #  J-3 

dX  -  dy2  -  A  ^ - 7 -  +  6 - T-  dx 


dW  -  dy2 


-  D  dyx 


dx2-  -  dy3 

cos  4> 


-  D'dy2 


D"  ,  2 

- -  dy 

cos  4 


-  gdx3 

3 

-  -*-dy3 

_i_3idy3 

cos  4  3X 


( IV-6U) 


The  metric  fundamental  tensor  is  given  therefore  by 


cosz4 


D'(D  +  — ^r)  D  %  + 

COS  4  cos  4 


D,{D  +  ^} 

cos  4 


D' 2  + 


cos24 


D'H  +  D-"-a 
1  cos  4 


(IV-65) 


D  H  +  P-  H?-- 
A  cos  4 


+ 

cos  4 


Hx  +  H2  +  7 


k.7  Malnardi-Codazzl  Equations;  Christoffel  Symbols 

It  is  now  very  easy  to  compute  Christoffel* s  symbols  of  the 
second  kind,  only  remembering  their  definition  as  given  in  (3.6): 

/  k  _h  3Yj  _X  3v^  _  3 vh  — 

I  ij  ^y1  3y  3  ay3  i 
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Before  doing  this,  we  must,  however)  assure  that  the  second 
derivatives  of  the  point  P  be  reversible,  i.e.  , 

a?i  /  3a  P  3*  P 

2y3  d  yi  *  dy^dy1  ^ 

This  leads  us  to  a  particular  form  of  the  so-called  Mainardi- 
Codazzi  equations  ( integrability  conditions)  for  the  coordinate 
surfaces. 


We  will  confine  oarsellres  to  the  computation  of  only  one  pair  of 
derivatives;  we  have  thus,  for  instance, 

8T’  ‘H*  '  '  <Jr  *  D,tg  *>4  -  -  D  sin  *  cos  *)  ♦  »•!,  , 


d\  3y* 


3-i  .  3Ji  .  . m  I  -,€2Ud»  tg»)—  ♦  D-i  , 

a*  ay1  d*  1  ^  'cos#  s» 


and,  therefore. 


♦  D  sin  #  cos  #  ♦  D"  tg  ♦  *  0  .  (IV-66) 

a*  ax 


Operating  in  a  similar  way  on  the  other  fundamental  vectors, 
we  easily  get  the  required  integrability  conditions: 
for  the  #  surfaces 


a  x  3i  * 


sin  ♦ 


0 


a  Ha  _  ay  i 
ax  aw  cos  * 


E,  sin  W  ♦  —  cos  ♦  ■  0 
g 


for  the  X  surfaces 

-  IS-  .  A  .  o 

a*  aw  g 

a  h,,  _  a  p'  i  , 

3*  aw  COS  t 


(IV-67) 
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for  the  W  surfaces 


3JL  -  iHL  ♦  D*  tg  *  -  0 

ax  3# 


3D"  aD 1 


♦  D  sin  9  cos  4>  +  Dn  tg  ♦  -  0  , 


3$  3  X 

The  Mainardi-Codazzi  equations,  or  integrability  conditions, 
simplify  in  tne  following  if  we  remember  the  meaning  of  and  Hg  : 

§r-&*v  t**-°  • 


3D" 

a  $ 


3D» 

ax 


♦  D  sin  $  cos  #  +  Dn 


tg  $ 


0 


3D  _  1  x  $  1 

aw  g  g  » 


( 17-68) 


3D1 

aw 


tg 


*3  1. 


3  1 

3^  g 


> 


ap"  _  i 

aw  "  g 


COS2*’  -  sin  $  cos  *  — 

3*  g 


t  JL 
“"ax7  «  * 


The  last  three  equations  yield  us  the  rate  of  variation  with 
the  dynamic  height  of  the  second  fundamental  tensor  of  the  equipotential 
surface,  and  allow  us,  therefore,  to  extend  in  the  third  dimension 
the  properties  of  our  equipotential  surface. 

It  should  be  noted  at  once  that  the  same  rate  of  variation  for 
the  intensity  of  gravity  is  yielded  by  Poisson-Brun' s  formula,  as  will 
be  seen. 


70 


WADC  TP.  52-149 


Moreover,  for  the  Christoffel  symbols  of  the  second  kind 


we  get: 


f  \x  ■  7  ( »■  ft  - w  ¥t  * D’  f*  * D  gh  cos  *)  > 

r\,  -A. '  'r  ( D"  It  ' D'  It  * D'  gh  cos  *)  • 

l  /  3D*'  3D'  3  D"  \ 

/-i  «  —  |  -  D,!  r —  +  2D"  -  D*  :r—  +  D"  gh  cos  #  )  -  sin  4>  cos  $  , 

/  A  V  3$  3X  2X  6  /  * 


^-l  i  l  /  3D  3D'  \ 

w* ** ) * 

r.,-r  32-KD"fT-D'lT)’  , 

fVH11" f^-D' It  C0S*-gJw  r 

r„ ■*(•”■«** lr-DlT*Dgk  cos*)’ 
T\,  ■  x  ( D  f; "  D’lr  * D’ gk  cos  *)  * tg  * 


3T)n  3D*  3D"  X 

D»  — - 2D'  —  +  D  “  +  D"  g  k  cos  *  )  , 

3*  3  X  <3X  / 


r.-K-ts-1' 

r-*  -r2  -  if  -  d* 

/23/33  »V  a*  Da»;> 


aH,  an,  i  .  .  , 

- *■  nnc  4.  _  Til  .=-=*■  -  S  "S’  *  k  COS  * 


(IV-69) 
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g  D;  r3  -  r*  *  -  g  d1;  r3  -  -  gD"  ; 
/  12  /  21  /  22 


rV/Vrvr  -0i  r 


g  222-*  (Hxh  +  Hgk  cos  *)  +  ~ 


g  (Dh3  +  2D*  hk  ♦  D«  k3)  + 


-  -2-  lg  g 
3W  S  K 


Also,  by  expanding  the  values  of  h  and  k  , 

d"^-d’ ■& > f- ■ 2 D,a ** * ]  > 

j"1  ■  7  [g  (  D»  -  »’  ^  Y  ♦  DD*  sta  *  coa  *]  , 

r1  -  7  fg  (  D"  ^  "  D'  -Jx  )  ~  +  T  (DD"  +  D'2)  sin  2  *  +  2  D"2  tg  #]  , 

I  32  8 

*(*"-***>  71  > 


7 [  (D"  ft  - s'  i  >  1 7  *  (■’”  *  *  4-  - 1 D' sin  2  *  4 )  t 


D’  cos  # 


-  j  U>"(A  -  *  a*  )  -  D'  (-fe  -  « 

■T  -D’  ^)f*2DD'tg*]  , 


(IT-70) 


r  -  T  tg  (D  &  -  D'  f  -  ?  sln  2  *  ]  ’ 
*  12 


D'  1 


WAHC  TR  52-149 


72 


.  i  [g  (D^L  -  D'^)  2L  -DD'sin  2  9  -  2  D'D"  tg  f  ]  , 


i  >»& 


D'  )  ■—  — 
39  d\  g 


(|  Dsin2*j#+D' 


+  —  COS2  9  ]  , 

g  * 


ilD(Jx 


3  g  v  ,  3  3  g  v  3g 

-  g  )-  O'  (J?  -  g  TT  )  , 


3*  '  3» 


-  g  D  i 


gD- 


-  g  D"  , 


iL  3*  g  -  2^  -  k*  f  *  _  .g. 
aw  g*  A 


lD"(^)a-2D-^|^+D(34)2] 


7a  • 

gA  cos*  9 


I*.  8  Components  of  Some  Vectors  and  Homographi.es 

It  is  now  easy  to  give  the  components  of  some  important  vectors 
and  homographies  in  our  fundamental  system  of  intrinsic  coordinates;  we 


get,  for  instance. 


for  the  vector  g  of  gravity 


v3  -  grad  W 


covariant  components:  *  W3 ■  §  3  *  (0,0,1)  , 

contravariant  components:  g*  ■  g.v*  -  T3.^  -  g3^  ■  (g31,  g32,  g33) 

,  g*h  cos  *  g2k  cos  9  _a\ 

“  .  »  “  .  »  g  t  9 


for  the  gravimetric  gradient  grad 
it  components:  gA  *  } 


covariant  components:  g/±  m 
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<V)|  O 


and  it  yields  us  the  rate  of  variation  of  the  intensity  of  gravity  with 
respect  to  the  dynamic  height. 


SECTION  V 


APPLICATIONS  OF  INTRINSIC  GEODESY 


5.1  Generalized  Theory  of  Botvos's  Torsion  Balance  (See  reference  53) 
We  consider  at  a  given  point  P  of  space,  referred  to  the  local 
astronomical  unit  vectors  (Ir),  a  beam  of  unit  length  symbolized  by  the 


unit  vector  «C: 


a  -  sin  V  cos  a  ix  *  sin  v  sin  a  i2  ♦  cos  v  i3  (V-l) 

(V  =  zenith  distance,  <i  *  azimuth  of  the  beam).  The  beam  may  be  sus¬ 
pended  at  P  at  its  baricenter;  the  acting  couple  will  therefore  be 
given  by  the  moment  vector  E  -  (&/a)  X  a,  2i/a  being  the  (vector) 
difference  between  the  gravity  at  the  ends  of  the  beam.  We  easily  get 
(Ua)X5  -  f-Tj-  (W  -  W22)  sin  2V  sin  a  +  T^g(sinaV  sin2  a-cos2v  )  + 

+  y  ( W13  sin2v  Sin  2  a  -  W1SJ  sin  2y  cos  a)|  Xi  + 

+  |i-(Wxl  -  W33)  sin  2  V  cos  o  +  W^cosV  -  sin2v  cos2  a)  + 

♦  y  (Wlg  sin  2v  sin2  a  -  W23  sin2V  sin  2  a)  j-  I2  +  (V-2) 

+  f—  (W_2  -  W  )  sin2V  sin  2  a  +  W  sin2V  cos  2  a  + 

1 2  22  n 

+  y  sin  2  V  (W23  cos  a  -  sin  a)  j  I3  , 

32  W 

where  W  stands  for  ,  and  the  xr  are  measured  along  the 

rs  3xr3xs 

axes  ir  . 
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We  thus  see  that  any  instrument  based  on  the  measure  of  moments, 
like  torsion  balances,  would  not  separate  the  second  derivatives  of  W 
with  equal  indices;  these  derivatives  may  only  be  determined  separately, 
remembering,  in  addition  to  the  results  of  the  balance  measurements, 
Laplace's  equation 


*  W33 


2  1112  (V-3) 

The  theory  of  all  types  of  balance  may  easily  be  derived  from 


the  foregoing  formulae  by  particularizing  the  direction  of  the  beam  and 
the  direction  along  which  the  scalar  moment  is  measured.  We  thus  get: 

EfttvBs  balance  of  the  first  kind  (horizontal  beam,  vertical 
swinging  axis)  , 


(tihr)  X  a  .  I3  -  1  (Waa  -  Wxi)  sin  2a  +  W1#  cos  2  a  ;  (V-fc) 


Eotvos  balance  of  the  second  kind  (oblique  beam,  vertical 
swinging  axis)  , 

ftj5)K5.I3  -  \  (W2a  -  Wxi)  sin*V  sin  2a  + 


(v-5) 


W12  sin2v  cos  2  a  +  -i-  sin  2  V  cos  a  -  Wi3  sin  a)  , 


Berroth's  balance  (oblique  beam,  horizontal  swinging  axis 
perpendicular  to  the  vertical  plane  containing  the  beam) : 

X  S.(-  sin  a  Tx  *  cos  a  ig)  «  7-  sin  2  V  [(  W13  -  W33)  cos3  a  ♦ 

(7-6) 

+  (W22  -W33)  sin2  a  +  Wis  sin  2a]  +  cos  2V(W13  cos  a  ♦  W33  sin  a  )  . 
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5.2 


We  already  have  generalized  in  1*.  2  on  an  arbitrary  surface 
Legendre's  developments  for  computation  of  geographical  coordinates 
and  azimuth.  We  are  now  able  to  further  generalize  the  same  develop¬ 
ments  in  space,  and  compute  both  the  geographic  coordinates  and  the 
dynamical  height,  assuming  that  the  path  is  a  straight  line  of  which 
the  azimuth,  the  zenith  distance,  and  the  length  are  given. 

We  immediately  have,  remembering  the  differential  equation  of 
a  straight  line  (3.8),  and  using  the  same  procedure  as  in  ii.2, 

(v-T) 


STJw  * « K  /*  -  -s M  ^ 


and  the  relationships  between  Xr  and  the  azimuth  and  zenith  distance 
of  the  path  are  given  by  the  formulae: 

X1  ■  -  —  sin  i  cos  a  +  ~  cos  9  sin  B  sin  a  +  —  (D’  ^  -  D"  — ^■■■)cos  z , 

A  A  A  oX  * 

( nr— 8) 

X*  ■  —  sin  z  cos  a  — 2.  COs  9  sin  z  sin  a  ♦  —  (D1  ^  -  D  6)  cos  z  , 

a  A  a  a#  ax 

X3  *  -  g  cos  z  , 


1  D'  X1  +  D"X*  ♦  X3  T 


cos  ♦ 


D  X1  +  D*  X2  +  X3  ^  — 

d9  g 


COS  B 


-i-  X3  . 

g 
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In  classical  Geodesy,  length  was  usually  measured  along 
geodesics  on  a  surface;  in  modern  Geodesy  there  is  a  tendency  to 
measure  length  along  geodesics  in  space  (straight  lines)  by  means  of 
electronic  or  interferometric  devices  (Shoran,  Hiran,  Geodimeter,  ...); 
the  foregoing  expansions  solve,  therefore,  the  first  problem  of  Geodesy 
in  this  case. 


Computation  of  Gravity  at  a  Point  of  Given  Intrinsic  Coordinates 


Let  us  have  two  points  in  space  PQ(y£)  and  P(yr),  and  put 

ayr  *  j*  -y£  (a/*  are  not  the  components  of  a  vector!);  let  us 
call,  moreover,  £  and  g^  the  gravity  vectors  at  P  and  PQ  respectively; 


we  then  will  have: 


g  “  gn  * 


*yi 


*  *  .  (V-9) 


Multiplying  the  foregoing  equation  by  the  fundamental  vectors 
vr  in  PQ,  we  get  the  covariant  components  gr  of  g  in  Pc  ;  we 


furthermore  remember  that 

3  w 

V'Tr-( 


ZJVi  •  vr  »  g 


i  •  vr  -  «r/i  “  «i/: 


,  _  r3 

i/r  *  ~  '  ri  * 


(V-10) 


Hi _ r 

9yi  dy-j 


vr  " 


.  V2 

6r/iJ  ‘  XT 


n\f  r  • 

'  rh  '  1  j  '  ih  '  r  j 


and  therefore 


r  r  l'riT  21  '  /rh'1J  /ih 'rH  (T_^j 

By  using  this  formula  we  may  compare  the  gravity  vectors  at  two 
different  points  of  space,  the  intrinsic  coordinates  of  which  are  given. 

If  only  the  value  g  of  the  intensity  of  gravity  is  required. 


we  simply  would  have 


g  + 

o 


»Ayi  +  21 


Ay1  +  •  •  •  (v-12) 
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5.J±  Computation  of  gravity  at  a  Point  of  Given  Polar  Coordinates. 

We  imagine  now  that  the  two  given  points  PQ  and  P  are 
connected  by  a  line  P  «  P(s)  of  length  sj  we  have  in  this  case 


8  -  1, 


♦  sim* 


si  l££i 

2  \ds7< 


♦  *  •  • 


(V-13) 


He  get 


32  -  _2£  dzL  -  3g  U 

ds  By1  ds  0yi  * 

45  -  ■?-£-,  xV  .  lH  x1/ 

ds*  3yi  0yj  dyi  i 


where  obviously  X*-  gives  us  the  components  of  the  tangent  unit  vectoi 
t  to  our  line  at  PQ,  and  depends  upon  the  nature  of  the  line. 

If  this  line  be  a  geodesic  (a  straight  line),  then  X^j  •  0  , 
and,  therefore#  we  have  for  a  straight  line 

-(Qi  •  S  (#-££ -r»  rXf  '*-■ 

In  a  perfectly  similar  way  we  get  for  the  intensity  of  gravity 


g  -  go 


*  •  •  •  •  (v-i5) 


In  all  of  the  above  formulae  the  X* ' s  are  related  to  the 
azimuth  and  the  zenith  distance  of  our  path  by  the  same  formulae  as 
in  5.2. 


The  Local  Cartesian  Equation  of  the  Geoid 
We  have  for  any  two  points  PQ  and  P 


p°  w 


bl-h 


a*  p 


Ar  3y% 


iy1  A  y«5 


(V-16) 
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We  remember  that 


SP 

V1 


'1  > 


B*  p 


(V-17) 


and,  therefore 

p  -  p0  -  |r(rj  ^r|AyiAy‘)  +  *  *  *  . 

By  multiplying  the  foregoing  formula  by  the  local  astronomic 
unit  vectors  (is  )  in  PQ  we  get 

(P-*o>-  is  -  K  *(?i-  UAyi  -  +  *  • 


(V-18) 


*, (V-19) 


where  the  scalar  expressions  of  v^*  ig  may  be  found  in  U.5j  we  thus 
have  the  cartesian  coordinates  referred  to  a  local  system  at  PQ  for 
a  point  of  given  intrinsic  coordinates.  If  we  take  Ay3*  0,  we  get 
the  cartesian  parametric  equation  of  the  level  surface  through  PQ . 

The  inverse  problem  of  getting  the  intrinsic  coordinates  from 
the  cartesian  has  been  already  solved  by  the  generalized  developments 
of  Legendre  given  in  5.2;  we  only  need  to  compute  previously  the  polar 
coordinates,  which  is  immediate. 


5.6  Application  to  Somigliana1  s  Field  (See  reference  hS,  U6,  50) 

It  is  well  known  from  Stoke1 s  theorem  that  a  gravitational 
field  is  fully  determined  in  the  empty  space  if  one  equipotential 
surface  of  the  field  and  the  total  acting  mass  (or  also  the  value  of 
gravity  in  one  point  of  the  surface)  are  given. 

Pizzetti  and  Somigliana  have  fully  solved  the  problem  of  de¬ 
termining  the  gravity  field  (gravitation  ♦  centrifugal  forces)  if  one 
of  the  equipotential  surfaces  is  an  ellipsoid  of  revolution.  Somigliana, 
in  particular,  has  furnished  a  formula  giving  the  value  of  gravity  along 
the  ellipsoid  itself. 
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The  formulae  of  Somigliana1  s  field  are  following  (on  the  ellipsoid) 


t(l  -  e2) 


,  N- 


(1-e2  sin2*)-^  "  (1-e*  sin2*) 


1 75 


a  ge  cos2*  +bgD  sin2* 
g  -  -p ===== _zr '  ,  (e  -  2  '» 

V  a2  cos2*  +  b2  sin2* 


(V-20) 


gp 


b.  e  +  i-  bw2  (1  +  e2  4  i2-  e4  ♦  * 

T  gs  2  w  v  35  180 


(Pizzetti*  s  function) 


)  . 


As  we  see,  only  the  knowledge  of  the  major  axis  a,  the  minor 
axis  b,  the  velocity  of  rotation  U).  and  the  gravity  ga  at  the  equator 
(equivalent  to  the  knowledge  of  the  total  mass)  are  required. 

The  formula  giving  the  value  of  g  has  been  expanded  into  a 
trigonometrical  series  of  the  latitude.  The  first  terms  of  the  series 
have  been  determined  by  Cassinis  and  3ilva,  and  give  the  international 
formula  for  gravity  accepted  by  the  International  Association  of  Geodesy 
at  the  Stockholm  meeting  in  1930. 

The  acceptance  of  an  ellipsoidal  field  as  a  standard  field  for 
gravity,  as  proposed  by  Somigliana  himself,  by  Dr.  ¥.  D.  Lambert,  and 
Professor  G.  Cassinis  has  eliminated  the  discrepancy  which  existed  up 
to  1930  between  geometrical  and  dynamical  Geodesy. 

The  international  formula,  based  on  the  International  Ellipsoid, 
is  the  following 

g  -  978, 0U9  (  1  ♦  0.005,288,4  sin2*  -  0.000,005,9  Sin2  2*)  . 

The  knowledge  of  gravity  along  the  ellipsoid  of  reference 
allows,  as  we  are  going  to  see,  the  easy  computation  of  all  elements 
needed  at  any  other  point  of  the  space. 
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We  have  on  the  basic  ellipsoid: 


^  -  P 


g  -  0 
la 


g 


13 


E  s 

g 


g 


ai 


g32  -  N2  cos2*  ,  g23  -  0  , 


(V-21) 


a31 


■£  .  a^g  e 


g  8  *  *  32 


8,«  -  0 


«33  * 


Moreover,  the  in tegr ability  conditions  we  have  found  give  us 
immediately 

8P._i.-jLi.  .  ^L.,1.  ♦tg*£..A  , 

o*  e  * 


aw 


g  a*2  g  *  3w 


g 


3N 


3^  "  <N  “  P)  tg  *  , 


( V-22) 


and  Laplace's  equation  gives,  furthermore. 


l«!  +  i  +  i_  ,aig_g  >a 

aw  g  g  >  n  ;  y,g  <  . 


(V-23) 


We  therefore  get  immediately  for  the  Christoff el  symbols  of  the 
second  kind: 

n-k*\  >  C  -  *  -  *£&■> , 

r' 


-iigjO. 

13  aw 

2  dig  N 


r 

1  23 


aw 


r1 «  a.  j^-ig-s  r 

/33  gp  a*  aw  *  '  1 


2 

12 


tg  *  , 


Fil-sp .  r„  -  e "  . 


(V-21») 


r. 


-  -  _gj£  ,fi 
33  aw 
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The  quantities  denoted  by  hs  k,  Hx  and  H#  have  the  following 


values : 


N  cos  ♦  31g  g 
g  3* 


.  .  i.  -?3&. S 


k  -  0  , 


H  -  0  . 


(V-25) 


g  3* 

,ts  of  th 

g/l  "  U'  "  “  ■§**  »  g/2  "  0  »  g/3  "  ”  g  rr 


The  covariant  components  of  the  gravimetric  gradient  are  given  by 

3  ,  (V-26) 

33 


and  the  fundamental  and  reciprocal  vectors  are 

»,  -»eo.*i;  ,rt  -i3  )  , 

*  -  7  I,)  .  **  - 

I1  3#  N  cos  # 


(V-27) 


The  foregoing  formulae  give  us  the  most  natural  generalization 
of  a  two-dimensional  geometry  on  the  ellipsoid,  as  usually  adopted  in 
Geodesy,  to  a  three  dimensional  scheme. 
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Colorado  Springs,  Colorado 

2  Commanding  General 
Strategic  Air  Command 

Offutt  Air  Force  Base,  Nebraska 

1  Commanding  General 

Tactical  Air  Command 
Langley  Air  Force  Base,  Virginia 

1  Commanding  General 

Technical  Training  Air  Force 
Gulfport.  Mississippi 

1  Commanding  General 

Air  Force  Missile  Test  Center 
Patrick  Air  Force  Base 
Cocoa,  Florida 

1  Commanding  General 

Air  Proving  Ground  Command 
Eglin  Air  Force  Base,  Florida 
ATTN:  Classified  Technical  Data 
Branch,  D/OI 

Commanding  General 
Air  Research  and  Dev.  Command 
P.0.  Box  1395 
Baltimore  1,  Maryland 

1  ATTN:  RDR 

1  ATTN:  RDK. 

1  ATTN:  RDT 

1  ATTN:  RDO 


Cyg  Activities 

5  Commanding  General 

Second  Air  Force 

Barksdale  Air  Force  Base,  La* 

2  Commanding  Officer 

USAF  Aeronautical  Chart  and 
Information  Service  Plant 
710  North  12th  Street 
St.  Louis  1,  Missouri 

2  Commanding  Officer 

USAF  Aeronautical  Chart 
and  Information  Service 
514  11th  St.,  N.W, 

Washington  25,  D.C. 

ATTN:  Technical  Library 

1  Commanding  General 

3415  Technical  Training  Wi^g 
Lowry  Air  Force  Base 
Denver,  Colorado 

1  Director 

Air  University  Library 

ATTN:  Req.  CR-3998 

Maxwell  Air  Force  Base,  Alabama 

5  Commanding  General 

Eighth  Air  Force 

Carswell  Air  Force  Base 
Fort  North,  Texas 

1  Air  Force  Eng.  Field  Representative 
Naval  Air  Missile  Test  Center 
Point  Mugu,  California 

1  Commanding  General 

Special  Weapons  Command 

Kirtland  Air  Force  Base,  New  Mexico 

1  Commanding  Officer 

Holloman  AF  Base 
ATTN:  6540th  Missile  Test  Wing 
New  Mexico 

1  Washington  AF  Eng.  Field  Office 

Room  4949,  Main  Navy  Bldg. 

Dept,  of  the  Navy 
Washington  25,  D.C. 


Cys  Activities 

1  Commanding  General 

Air  Force  Cambridge  Res.  Center 
230  Albany  Street 
Cambridge  39,  Massachusetts 

5  Commanding  General 

Fifteenth  Air  Force 
March  Air  Force  Base 
California 

3  Commanding  General 

5th  Air  Division 
APO  118,  C/0  Postmaster 
New  York,  New  York 

3  Commanding  General 

7th  Air  Division 
APO  125,  C/O  Postmaster 
New  ^ork.  New  York 

Army 

1  Chief  of  Army  Field  Forces 

Fort  Monroe,  Virginia 

1  Commanding  Officer 

White  Sands  Proving  Ground 

Oro  Grande,  New  Mexico 

1  Commandant 

Army  War  College 

Fort  Leavenworth,  Kansas 

1  Chief  of  Engineers 

ATTN:  Res.  and  Dev.  Division 
Building  T-7,  Gravelly  Point 
Washington  25,  D.C. 

1  Commanding  Officer 

Engineer  Res.  and  Dev.  Labs. 

The  Engineer  Center 
Fort  Belvoir,  Virginia 

1  Commanding  Officer 

Army  Map  Service  Library 
Corps  of  Engineers 
6500  Brooks  Lane 
Washington  16,  D.C. 

1  Commandant 

National  War  College 
Washington  25,  D.C. 
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C£8  Activities 
Navy 

1  Director 

Special  Devices  Center 
Office  of  Naval  Research 
ATTN:  Visual  Design  Branch  940 
Sands  Point 

Port  Washington  L.I.,  N,Y. 

1  Hydrographer 

U„S0  Navy  Hydrographic  Office 
Department  of  the  Navy 
Washington  25,  D.C. 

1  Officer  in  Charge 

U.S.  Naval  Photographic 
Interpretation  Center 
Naval  Receiving  Station 
Washington  25 ,  D.C, 

Other  U.S,  Government  Agencies 

1  Central  Intelligence  Agency 

ATTN:  Office  of  Collector  of 

Documentation,  Control 
No„  CD-A-18831 
2430  1  Street,  N.W. 

Washington  25,  D.C. 

1  Director 

U.S,  Coast  and  Geodetic  Survey 
Commerce  Department 
Washington  25,  D.C, 


Cys  Activities 


1  U.S,  Geological  Survey 
1033  Interior  Bldg. ,  N. 

Washington  25,  D.C, 

Others 

10  Ohio  State  University  Research  Pound. 
Mapping  and  Charting  Research  Lab. 
2593  West  Hardin  Road 
Columbus  10,  Ohio 

1  North  American  Aviation,  Inc. 

ATTN:  Aerophysic8  Library 
12214  Lakewood  Blvd. 

Downey,  California 

1  Northrop  Aircraft,  Inc. 

ATTN:  Mr.  John  Northrop 
Hawthorne,  California 

1  RAND  Corporation 
1500  4th  Street 
Santa  Monica,  California 
Thru:  WCRR 
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